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SECTION  1 


INTRODUCTION 


Silicon  microcircuits  which  use  Si02  for  the  gate  dielectric  and/or  surface 
passivation  degrade  when  exposed  to  ionizing  radiation.  The  degradation  results  from 
two  phenomena: 

1.  Radiation  induced  trapped  holes,  A  Not',  and 

2.  Radiation  induced  interface  states  A  Nit- 

The  details  of  the  creation  and  annihilation  of  these  two  charge  components  have 
been  studied  for  the  past  30  years  since  the  discovery  of  ionization  induced  degrada¬ 
tion  of  transistors  in  an  early  satellite.  The  research  has  involved  hole-electron  pair 
generation  and  recombination  in  Si02,  hole  transport  and  trapping,  interface  state 
formation  mechanisms,  the  chemical  structure  of  the  hole  trap  and  interface  state,  the 
energy  density  of  the  interface  states,  the  detrapping  mechanisms  of  trapped  holes, 
and  the  annihilation  mechanisms  f  interface  states.  Many  models  have  been  devel¬ 
oped  for  the  observed  characteristics  and  their  dependence  on  dose,  dose  rate,  electric 
field,  temperature,  oxide  thickness,  and  oxide  processing  details.  Most  of  the  re¬ 
search  to  characterize  the  radiation  induced  trapped  hole  and  interface  state  densities 
has  been  done  with  a  nearly  ideal,  one-dimensional  test  structure,  the  MOS  capacitor, 
MOSC.  The  primary  techniques  used  to  study  the  mechanisms  of  A  Not  and  A  Njt  on 
capacitors  involve  capacitance-voltage  or  conductance-voltage  characteristics.  These 
measurements  are  made  at  various  temperatures,  sweep  rates,  and  frequencies  to  char¬ 
acterize  the  trapping  and  detrapping  of  ANot  and  the  energy  dependence  of  ANit . 
Other  techniques  involving  charge  injection  and  gate  current  measurements  are  used 
to  measure  the  spatial  dependence  and  energy  levels  of  the  trapped  holes. 

Although  the  MOSC  is  an  ideal  structure  for  studying  the  physics  of  ANot 
and  ANit,  it  is  not  a  practical  test  structure  for  engineering  applications  such  as 
device  characterization  and  hardness  assurance.  Some  of  the  reasons  include: 
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1.  MOSCs  are  not  used  to  build  microcircuits. 


2.  MOSCs  cannot  be  used  to  study  parasitics,  e.g.,  sidewall  and  backchannel  leak¬ 
age  in  SOS  and  SOI. 

3.  Practical  bias  conditions  cannot  be  evaluated,  e.g.,  PN  junction  bias. 

4.  Geometrical  effects  cannot  be  evaluated,  e.g.,  gate  length  and  width. 

5.  MOSCs  are  usually  large  and  hence  are  seldom  included  on  microcircuit  die  as 
test  structures. 

It  has  long  been  recognized  that  the  most  practical  test  structure  for  characterizing 
the  ionization  response  of  MOS  microcircuits  is  the  MOS  transistor,  MOSFET,  since 
it  is  the  primary  circuit  element  which  determines  the  failure  level  of  the  microcircuit. 
It  has  also  been  shown  that  in  bipolar  microcircuit  technologies  which  use  recessed 
field  oxides,  the  dominate  failure  mode  (buried  layer  to  buried  layer  leakage)  is  best 
characterized  by  a  parasitic  MOSFET  structure.  There  are  two  principal  reasons  why 
the  MOSFET  has  not  been  used  extensively  in  the  past  as  the  test  structure  of  choice 
to  study  and  characterize  A  Not  and  A  Nit-  First,  the  techniques  for  separating 
the  two  charge  components  using  transistor  measurements  have  only  recently  been 
developed.  Second,  the  characteristics  which  make  the  transistor  the  more  practical 
test  structure  often  result  in  limitations  which  make  it  less  ideal  for  studying  the 
physics  of  ANot  and  AJV/r.  For  example,  the  geometric,  bias  and  parasitic  effects 
which  are  often  important  in  understanding  failure  mechanisms  in  actual  microcircuits 
often  result  in  a  distortion  of  the  transistor  electrical  characteristics  which  make  it 
difficult  to  separate  out  the  A  Not  and  A  Nit  in  the  active  gate  oxide. 

If  the  MOSC  is  the  ideal  structure  for  studying  the  physics  of  ANot  and 
A  Nit,  and  the  MOSFET  is  the  ideal  structure  for  studying  the  component  degrada¬ 
tion  which  results  in  microcircuit  failure,  then  it  is  reasonable  to  question  why  charge 
separation  needs  to  be  performed  on  MOSFETs. 

Since  the  development  of  MOS  microcircuit  technologies  in  the  late  1960’s, 
most  analyses  of  microcircuit  failure  were  based  on  the  threshold  voltage  shift,  AVth, 
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which,  in  early  MOS  technologies,  had  a  negative  shift  due  to  trapped  holes.  How¬ 
ever,  as  process  hardening  techniques  developed  to  reduce  the  hole  trapping  efficiency 
and  oxides  were  thinned,  resulting  in  fewer  available  holes  per  rad,  the  trapped  hole 
density  was  no  longer  the  only  charge  component  which  had  to  be  considered.  Several 
failure  mechanisms  relating  to  interface  states  started  to  become  dominate  such  as, 
a)  mobility  degradation,  and  b)  positive  threshold  voltage  shifts  in  n- channel  MOS- 
FETs  referred  to  as  “rebound”  or  “super-recovery”. 1  Thus,  in  many  modern  MOS 
microcircuit  technologies  which  have  thin  oxides,  or  have  been  process  hardened  to 
reduced  hole  trapping,  both  A  Not  and  AN  it  are  important.  Also,  because  A  Vor>  the 
voltage  shift  resulting  from  A  Not*  and  AV/r,  the  voltage  shift  resulting  from  A  Nit, 
have  opposite  signs  in  n-channel  devices,  one  can  perform  “hardening”  by  balancing 
the  two  components.  However,  since  the  buildup  and  decay  of  AVot  and  AV/r  have 
different  time  constants  associated  with  them,  the  balancing  of  AVor  and  AV/r  only 
works  over  a  limited  range  of  dose  rates  and/or  annealing  times.  Also,  buildup  of  in¬ 
terface  states  can  degrade  mobility  resulting  in  speed  degradation  of  the  microcircuit. 
Therefore,  “total  dose”  response  is  not  single  valued,  since  the  valu^i  of  AVqt  and 
AV/r  at  a  given  total  dose  may  vary  significantly  depending  on  the  dose  rate  or  mea¬ 
surement  time.2  It  is  mandatory  in  most  modern  MOS  technologies  to  separate  AVor 
and  AV/r  in  order  to  evaluate  the  ionization  hardness.  A  part  which  is  hard  to  many 
megarads  when  tested  with  a  60Co  source  at  200  rads/s,  one  hour  after  irradiation 
(see  MIL-STD-883  method  1019),  may  fail  at  a  significantly  lower  dose  when  tested 
at  high  dose  rate  (weapon  environment)  or  low  dose  rate  (space  environment).  The 
reason  is  that  the  part  may  have  a  high  AVor  and  a  high  AV/r  at  intermediate  dose 
rates  which  balance  to  prevent  failure.  At  the  high  dose  rate,  failure  is  dominated  by 
A  Not  mechanisms,  and  at  low  dose  rates  failure  is  dominated  by  A  Nit  mechanisms. 
In  fact,  the  values  of  AVor  and  AV/r  could  possibly  be  used  as  hardness  measures 
for  the  high  dose  rate  and  low  dose  rate  environments,  respectively. 

Thus,  there  are  many  practical  engineering  applications  for  charge  separa¬ 
tion  in  MOSFETs,  among  them: 

1.  Characterization  of  the  hardness  of  a  technology. 
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2.  Evaluation  of  the  effectiveness  of  process  techniques  to  reduce  ANot  and  A  Nit- 

3.  Monitoring  the  process  to  control  ANot  and  A  Nit- 

4.  Qualification  and  hardness  assurance  of  process  lots. 

For  engineering  applications,  such  as  the  four  mentioned  here,  it  is  probably 
not  necessary  to  determine  the  actual  values  of  ANot  o r  ANjt,  but  rather  the  voltage 
components  AVot  and  AV1T.  A  For  is  proportion'd  to  the  “effective”  trapped  hole 
density  and  AVIT  is  a  measure  of  the  net  charge  of  the  interface  states  at  threshold. 

There  have  been  a  number  of  charge  extraction  techniques  developed  for 
MOSFETs  over  the  past  few  years.  However,  most  of  these  techniques  have  been 
developed  to  extract  ANjt(AV!t).  Only  one  technique  has  been  proposed  to  extract 
ANot{AVot).3  The  ANot  technique  is  based  on  the  assumption  that  the  interface 
states  are  predominantly  acceptors  in  the  upper  half  of  the  Si02  band  gap  and  donors 
in  the  lower  half  and  that  the  net  charge  of  the  interface  states  is  approximately 
zero  at  midgap.4  This  assumption  has  long  been  used  to  extract  ANIT  using  MOSC 
C-V  characteristics.  Several  techniques  have  been  developed  to  extract  A NIT  from 
MOSFET  characteristics.  Among  these  are: 

a.  Subthreshold  swing;6,6 

b.  Weak  inversion;7 

c.  Mobility;8 

d.  Charge  pumping;9  and 

e.  Deep  level  transient  spectroscopy.10 

The  weak  inversion  technique  is  used  to  find  the  value  of  A  Nit  at  a  single 
surface  potential  half  way  between  midgap  and  inversion.  The  mobility  technique  is 
based  on  the  work  of  Sun  and  Plummer11  who  showed  that  the  process  induced  oxide 
fixed  charge,  Qss,  which  is  located  within  50  A  of  the  Si02— Si  interface,  caused  a 
degradation  of  surface  mobility  according  to  the  relation 

_  Mo 

^  1  +  ocQss  * 
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where  n  is  the  degraded  mobility,  fi0  the  “ideal”  mobility,  and  a  an  empirical  coeffi¬ 
cient.  Galloway8  used  this  relationship  to  describe  the  radiation  induced  degradation 
of  mobility  by  interface  states  in  the  form 

=  Mo 

^  1  +  ttAJV/y  * 

where  /z  is  now  the  post-irradiation  mobility,  n0  the  pre-irradiation  mobility,  and 
a  is  an  empirical  coefficient  which  may  vary  with  process.  Charge  pumping  is  a 
technique  which  can  be  used  to  measure  A Njt  based  on  a  substrate  (channel)  current 
proportional  to  the  number  of  interface  states  charged  and  discharged  by  sweeping  the 
surface  potential.  The  charge  pumping  technique  requires  a  separate  contact  to  the 
substrate/channel  and  does  not  use  conventional  MOSFET  I-V  characteristics.  The 
deep  level  transient  spectroscopy  (DLTS)  technique  has  only  recently  been  applied 
to  MOSFETs.  This  technique  is  not  practical  for  engineering  applications  because 
measurements  must  be  made  over  a  wide  range  of  temperatures. 

The  subthreshold  technique,  which  was  first  presented  in  1984s  and  was  for¬ 
malized  in  1986, 12  is  based  on  standard  I-V  characteristics  and  is  the  only  technique 
which  combines  a  method  for  AVbr  and  AV/j*.  The  subthreshold  drain  current  at  a 
fixed  drain  to  source  voltage,  Vos,  is  measured  as  a  function  of  gate  voltage  from  the 
leakage  current  (or  limiting  resolution  of  the  measurement  system)  through  inversion. 
In  an  ideal  device,  the  drain  current  and  gate  voltage  are  related  by  Id  oc  exp  (Vb). 
When  plotted  as  log  Id  versus  Vb,  the  straight  I-V  characteristic  can  be  extrapolated 
to  a  calculated  midgap  current.  Comparing  the  pre-  and  post-irradiation  character¬ 
istics,  the  midgap  voltage  shift,  AVmg,  as  well  as  the  change  in  subthreshold  swing 
(inverse  slope),  AS,  can  be  determined.  The  value  of  A Vmq  is  equivalent  to  AVbr 
and  AS  is  proportional  to  A Vrr.  The  subthreshold  charge  separation  technique  has 
proven  to  be  the  easiest  to  perform  and  is  the  most  widely  used.  However,  as  will  be 
discussed  in  detail  in  this  document,  there  are  many  limitations  to  the  technique. 

The  purpose  of  this  document  is  to  provide  guidance  in  the  use  of  the  sub¬ 
threshold  charge/voltage  separation  technique.  To  this  end,  the  theory  and  assump¬ 
tions  are  discussed  in  detail.  The  applications  of  the  technique  are  discussed  with  re¬ 
spect  to  electrical  measurement  techniques  and  data  analysis.  Finally,  the  limitations 
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of  the  technique  are  discussed  with  regard  to  interferences  from  leakage,  parasitics, 
and  violation  of  assumptions.  In  the  Appendix,  a  proposed  measurement  and  data 
analysis  standard  is  given. 
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SECTION  2 


MOSFET  SUBTHRESHOLD  THEORY 


MOSFET  subthreshold  theory  refers  to  the  electrical  characteristics  of  the 
transistor  in  weak  inversion.  In  weak  inversion,  the  current  in  the  channel  is  dom¬ 
inated  by  diffusion  of  minority  carriers.  This  differs  from  strong  inversion  channel 
current  which  is  dominated  by  electric  field  assisted  majority  carrier  transport.  The 
subthreshold  current  is  used  to  determine  the  shift  in  threshold  voltage  due  to  trapped 
holes,  AVot,  and  interface  states,  A V/r*  The  accuracy  of  this  calculation  is  depen¬ 
dent  on  whether  assumptions  used  to  develop  the  subthreshold  charge  separation 
technique  are  violated.  The  theory,  together  with  the  assumptions,  predicts  a  near 
linear  subthreshold  swing  (mV/decade)  for  a  MOSFET.  A  nonlinear  swing  may  imply 
a  violation  of  an  assumption  or  the  effect  of  a  variable  not  considered  in  the  theoreti¬ 
cal  development.  In  this  section,  the  MOSFET  subthreshold  theory  is  presented.  In 
Section  3,  violations  of  assumptions  are  checked  to  determine  their  theoretical  effect 
on  the  subthreshold  charge  separation  technique.  Most  important,  their  effect  on  the 
subthreshold  swing  and  analysis. 


2.1  THEORETICAL  MOSFET  WEAK  INVERSION  CHANNEL 
CURRENT. 

In  a  MOSFET,  the  theoretical  subthreshold  channel  current,  /0,  is  given  by 

where 


H  =  effective  mobility, 

N  =  interface  doping, 

n,-  =  intrinsic  doping  concentration, 
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Z  =  channel  width, 
L  =  channel  length, 


fi  =  q/kT, 

(j>s  =  surface  potential, 
a  =  y/2  ( s */ Ld)  /C0*, 


e,  =  silicon  dielectric  constant, 
Ld  =  debye  length, 

Cox  =  toxl^oxt 

eox  =  oxide  dielectric  constant, 
foa:  =  oxide  thickness. 


and  drain  voltage  is  greater  than  a  few  kT/q. 


This  equation  is  derived  using  Brews’  charge  sheet  model.13,6  The  current 
in  this  case  refers  to  the  current  which  flows  in  the  channel  of  the  device.  It  does  not 
include  gate  and  substrate  leakages  or  any  other  parasitic  effect.  The  channel  current 
in  equation  (1)  is  a  strong  function  of  the  surface  potential.  Gate  voltage  is  related 
to  the  surface  potential  so  that  the  dVa/dinlu  can  be  calculated.  This  derivative 
gives  the  gate-voltage  swing  needed  to  reduce  the  current  one  decade,  and  is  called 
the  subthreshold  swing.  The  subthreshold  swing,  S,  is  defined  by 


S  =  p-HnlO 


dpVg  dp<f>S 
dfi4>s  dlnl0 ' 


(2) 


Changes  in  the  variables,  which  determine  the  channel  current,  can  affect 
the  subthreshold  swing,  if  they  are  a  function  of  the  surface  potential.  Equation  (l) 
is  used  to  determine  the  midgap  current  in  the  subthreshold  technique.  It  can  also 
be  used  to  determine  the  inversion  current  and  thereby  the  inversion  voltage,  if  the 
threshold  voltage  is  not  measured.  Changes  in  any  of  the  variables  in  equation  (1) 
will  change  the  value  of  the  midgap  and  inversion  currents. 

To  calculate  the  derivative  in  equation  (2),  a  function  relating  the  gate  volt¬ 
age  to  the  surface  potential  is  needed.  Using  depletion  charge,  the  following  equation 
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was  constructed  by  Brews,13 


0Va=Ms+a(Ms-  1)1/S- 


(3) 


Using  equations  (2)  and  (3),  and  assuming  all  coefficients  are  constant  with 
surface  potential,  the  derivative  in  equation  (2)  can  be  calculated  and  the  subthreshold 
swing  becomes 


S  =  p-'lnlO 


(4) 


where 

CD  —  channel  depletion  capacitance. 


This  is  the  equation  on  which  the  subthreshold  charge  separation  technique 
is  based.  Usually  the  denominator  in  equation  (4)  is  assumed  to  be  1.0  (i.e.,  it  is 
assumed  that  a  »  CD/C0X). 


The  subthreshold  swing  is  not  theoretically  constant  from  midgap  to  inver¬ 
sion.  This  is  due  to  the  depletion  capacitance  terms  in  equation  (4).  The  depletion 
capacitance  results  from  the  depletion  of  majority  carriers  in  the  channel  region  and 
is  in  series  with  the  insulator  (gate)  capacitance.  As  the  surface  potential  is  increased, 
the  depletion  region  grows  and  the  total  capacitance  decreases.  The  total  capacitance 
continues  to  decrease  until  an  inversion  layer  is  formed  and  the  total  capacitance 
sharply  increases.  This  process  affects  the  subthreshold  swing,  because  the  depletion 
capacitance  varies  with  the  surface  potential.  The  minimum  total  capacitance  occurs 
between  midgap  and  inversion,  and  this  is  also  the  point  at  which  the  minimum  sub¬ 
threshold  swing  occurs.  The  swing  increases  as  the  surface  potential  is  increased  or 
decreased.  For  a  typical  device,  the  pre-irradiation  deviation  from  linear  is  small,  as 
is  shown  in  Figure  1. 
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Figure  1.  Difference  of  the  theoretical  subthreshold  drain  current  from 
linear  due  to  depletion  capacitance. 


2.2  SUBTHRESHOLD  CHARGE  SEPARATION  TECHNIQUE. 


There  are  two  methods  to  separate  the  threshold  or  inversion  voltage  into 
a  voltage  shift  due  to  oxide  trapped  holes,  AV$r>  and  interface  states,  A Vjy,  using 
the  subthreshold  technique.  The  two  methods  depend  on  whether  AV0r  or  A Vjt  is 
directly  calculated  from  the  subthreshold  curve  and  are  referred  to  as: 

1.  Midgap  voltage  method, 

2.  Subthreshold  swing  method. 

The  application  of  either  method  is  dependent  on  whether  the  measured 
threshold  voltage  or  the  calculated  inversion  voltage  is  used. 

2.2.1  Determination  of  AVinv  and  A Vw 

2 .2 . 1.1  Inversion  Voltage.  The  onset  of  strong  inversion  occurs  at  a  gate  voltage 
corresponding  to  a  surface  potential,  <j>s,  equal  to  twice  the  bulk  fermi  level,  <f>B,  where 

4,  =  Hb  =  —In  (-)  .  (5) 

The  gate  voltage  at  this  surface  potential  is  referred  to  as  the  inversion  voltage,  VjyV. 
The  channel  current  for  strong  inversion,  Imv,  is  calculated  using  equation  (1)  with 
<f>s  =  2 4>b-  Vjnv  is  the  gate  voltage  at  Ijyv  which  is  determined  by  interpolating 
the  Id  versus  Vq  data.  Ijsv  and  Vjsv  are  illustrated  in  Figure  2.  Figure  2  is  a  plot 
of  log  Id  versus  VG  at  VD  =  100  mV  for  a  typical  n-channel  MOSFET.  In  order  to 
calculate  Imv »  the  following  parameters  must  be  known:  Z,  Leff,n,tox,  N,  and  Vp. 
This  parameter  set  can  be  reduced  if  the  peak  transconductance,  Gm ,  is  determined 
from  the  experimental  data.  Since 


Gm  = 


AId 
A Vq 


Vd 


ZnCozVD 

L.ff  ’ 


(6) 
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then 


Vp  >  kT/q , 


(7) 


/  =  Gm  a  n,-2  ew* 

""  Vd  2/?2  iV  (20<j>B  -  if3 

where  Vp  is  the  drain  voltage. 

Using  equation  (7),  the  parameter  set  required  to  calculate  ijjvv  is  reduced 
to  tox  and  N,  since  the  remaining  information  can  be  determined  from  the  experimental 
data. 


The  pre-irradiation  value  of  the  inversion  voltage,  Vjsvu  is  found  from  the 
pre-irradiation  value  of  Iinv  using  the  pre-irradiation  value  of  Gm-  However,  because 
of  the  exponential  dependence  of  Id  on  VG,  a  factor  of  2  or  3  error  in  IiNv  will  only 
result  in  a  small  percentage  change  in  Vinvi •  The  error  will  remain  small  if  the  post¬ 
irradiation  value  of  IINV  is  set  equal  to  the  pre-irradiation  value.  Therefore,  the  post¬ 
irradiation  value  of  the  inversion  voltage,  ViNy^  can  be  determined  by  interpolating 
the  post-irradiation  ID  versus  VG  curve.  The  irradiation  induced  inversion  voltage 
shift,  AVinv,  is  determined  from  the  following  relation 


A Vinv  =  Vinvi  ~  Visvb  (8) 

2. 2. 1.2  Threshold  Voltage.  The  parameter  most  often  used  to  characterize  the 
radiation  hardness  of  an  MOS  transistor  is  the  threshold  voltage  VTfi,  For  drain 
voltage  greater  than  {yG  —  VTH ),  the  transistor  is  in  saturation,  and  Ip  oc  (VG  —  Vth)2- 
Therefore,  a  plot  of  versus  VG  should  yield  a  straight  line.  If  \JTp  versus  VG  is 
extrapolated  to  ID  =  0,  the  VG  intercept  will  be  Vjh •  If  the  drain  voltage  is  less  than 
[VG  —  Vj7f),  the  transistor  is  in  the  linear  region,  and  a  plot  of  Ip  versus  VG  should 
yie'd  a  straight  line.  VTH  is  found  by  extrapolating  Ip  versus  VG  to  Ip  =  0.  AVTh  is 
given  by 

AVTh  =  VTHy  -  VT:tI,  p) 

where  Vthi  and  Vthi  are  the  post-  and  pre-irradiation  threshold  voltages,  respectively. 
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2.2.2  Midgap  Voltage  Method. 


The  midgap  voltage  method  of  the  subthreshold  charge  separation  technique 
is  based  on  two  calculations:  l)  the  radiation  induced  voltage  shift  from  oxide  trapped 
charge,  AV0t,  using  the  change  in  midgap  voltage,  A Vmg,  and  2)  the  radiation  in¬ 
duced  voltage  shift  from  interface  states,  A V/r,  using  A Vmg  and  the  change  in  the 
threshold  voltage,  A Vth,  or  the  inversion  voltage,  AV/jw. 


2. 2.2.1  A  Vqt-  The  value  of  AVbr  is  obtained  from  A  Vmg  assuming  the  following 
relation: 

AVer  =  A  VMc  =  (10) 

where  Cox  —  Eoz/tox  =  3.45  x  10 19/tox  F/cm  for  Si02,  and  A  Not  is  the  radiation 
induced  trapped  hole  density.  Equation  (10)  is  based  on  the  following  assumptions: 


a.  The  trapped  hole  density  is  a  uniform  sheet  charge  at  the  interface;  and 

b.  The  net  charge  of  the  interface  states  at  midgap  is  zero.4 


For  relatively  thin  gate  oxides  (200-1000  A),  the  spatial  distribution  of  the 
trapped  holes  has  been  determined  using  a  number  of  indirect  techniques.  Except 
for  a  narrow  region  within  20-50  A  of  the  interface  where  trapped  holes  are  quickly 
compensated  or  removed  by  electrons  tunneling  from  the  silicon,  the  trapped  holes  are 
located  within  50-150  A  of  the  oxide.14  The  actual  spatial  dependence  of  the  trapped 
hole  density  within  this  region  is  not  well  known  but  is  probably  best  represented  by 
a  double  exponential  peaking  around  30-50  A.  For  some  oxides,  e.g.,  commercial  soft 
gate  oxides14  and  thick  field  oxides,15  the  holes  may  be  nearly  uniformly  trapped  in 
the  oxide  bulk.  In  the  case  of  the  uniform  trapping,  the  midgap  shift  is  only  1/2  of 
the  value  for  a  sheet  charge  at  the  interface  and  hence, 


avmg  = 


qANor 

~2 C^T' 


(11) 
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In  the  application  of  the  subthreshold  charge  separation  technique,  the  ra¬ 
diation  induced  interface  states  are  assumed  to  be  primarily  acceptors  in  the  upper 
half  of  the  band  gap  and  donors  in  the  lower  half.  Near  midgap,  the  net  charge  of 
the  interface  states  is  approximately  zero.4  If  the  net  charge  of  the  interface  states  at 
midgap  is  zero,  then  the  change  in  midgap  voltage  is  a  result  of  the  net  oxide  trapped 
charge  which  is  assumed  to  be  holes. 

The  midgap  voltage  usually  cannot  be  measured  directly  because  the  chan¬ 
nel  current  at  midgap  is  usually  smaller  than  the  leakage  current  and/or  too  small 
to  measure  with  conventional  picoammeters  or  electrometers.  Therefore,  an  extrap¬ 
olation  technique  is  used  to  determine  the  midgap  voltage  from  channel  current, 
versus  gate  voltage,  VG,  in  sub  threshold.  The  channel  current  in  sub  threshold  (weak 
inversion)  is  given  in  equation  (l). 

The  channel  current  at  midgap,  Img ,  can  be  found  by  substituting  the 
midgap  value  of  surface  potential,  4b,  where 

ts  =  *s  =  —In  (-)  ,  (12) 

q  \riiJ 

into  equation  (1).  The  midgap  voltage,  VMG  is  the  gate  voltage  at  Img  found  by 
extrapolating  the  measured  subthreshold  characteristic.  Figure  2  is  a  plot  of  loj:  Id 
versus  VG  at  VG  =  100  mV  for  a  typical  n-channel  MOSFET  showing  the  extrapolation 
of  Id  to  Img-  Using  equation  (1),  Img  is  calculated  similar  to  J/jvy  with  <j>s  =  <f>B- 

The  pre-irradiation  value  of  midgap  voltage,  Vmgi,  is  found  from  the  pre¬ 
irradiation  value  of  Img  using  the  pre-irradiation  value  of  Gm-  The  value  of  Img 
calculated  in  this  manner  will  be  somewhat  low,  since  the  mobility  in  subthreshold 
may  not  be  constant  and  may  increase  as  the  surface  potential  is  lowered.  Therefore, 
the  post-irradiation  value  of  midgap  voltage,  Vmg-v  can  be  determined  by  extrapolat¬ 
ing  the  post-irradiation  Id  versus  VG  curve  to  Imgi-  A Vmg  is  determined  from  the 
following  relation 

A  VMg  =  Vmg-i  ~  Vmgi-  (13) 
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2. 2. 2.2  A Vit-  The  radiation  induced  voltage  shift  due  to  charged  interface  states 
is  assumed  to  be  that  portion  of  the  inversion  voltage  shift  or  threshold  voltage  shift 
which  is  not  due  to  the  trapped  holes.  The  use  of  inversion  voltage  is  more  mathe¬ 
matically  correct  since  it  is  the  value  of  gate  voltage  for  a  specified  value  of  surface 
potential.  However,  the  use  of  threshold  voltage,  which  is  an  empirical  value,  is  more 
practical  since  it  is  the  parameter  specified  on  data  sheets  and  generally  used  to  mea¬ 
sure  radiation  performance. 

The  inversion  voltage  change  due  to  irradiation,  AV/^y,  is  a  result  of  the 
voltage  shift  due  to  oxide  trapped  charge,  A Vor>  and  the  voltage  shift  due  to  the 
charged  interface  states,  A  V/r-  In  n-channel  devices,  the  value  of  A  Vit  is  positive  due 
to  the  predominance  of  acceptors,  and  in  p-channel  devices,  A  Vjt  is  negative  due  to 
the  predominance  of  donors.  Using  the  midgap  voltage  method,  the  value  of  A  Vit  is 
found  from  the  following  relation: 


AV/t  =  AVW  -  AVor.  (14) 

As  with  the  midgap  voltage  shift,  the  pre-irradiation  and  post-irradiation  values  of 
Ijnv  are  calculated  using  the  pre-irradiation  value  of  transconductance. 

Using  AVTh  instead  of  AV/^y,  A V/y  is  determined  from  the  following  rela¬ 
tion: 

A  Vjt  —  AVjt/  —  A  VqT'  (15) 


2.2.3  Subthreshold  Swing  Method. 

The  subthreshold  swing  method  of  the  charge  separation  technique  is  based 
on  two  calculations:  1)  the  radiation  induced  voltage  shift  due  to  interface  states, 
A  Vjt,  using  the  change  subthreshold  swing,  AS,  and  2)  the  radiation  induced  ohift 
from  oxide  trapped  holes,  AVor>  using  A  Vjt,  and  the  change  in  threshold  voltage, 
AVjh,  or  the  inversion  voltage,  A Vinv 
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2. 2.3.1  A Vjt.  Assuming  a  »  C/j/Cox,  equation  (4)  can  be  rewritten,  which  gives 


S  = 


kT 

— £nlO(l  +  CD/C0X) . 


(16) 


If  the  pre-irradiation  value  of  the  subthreshold  swing,  S/,  is  defined  by 
equation  (16),  the  post-irradiation  subthreshold  swing,  is  defined  by 


kT 

S~  =  —InlO 

q 


i  + 


Cd  +  C/r 


(17) 


The  difference  between  the  pre-  and  post-irradiation  rubthreshold  swings,  AS,  is 


kT  C1T 

A  S  =  Sn-St=  — lnlO—p~. 

Q  C„x 


(18) 


If  the  interface  states  are  uniform  between  midgap  and  inversion,  the  interface  state 
capacitance,  C/r,  causes  the  voltage  shift,  AV/r,  which  is  equivalent  to 


A  Vjt  =  GitIC0X<I>b  (19) 

Using  equations  (18)  and  (19),  the  value  of  A  Vjt  is  obtained  from  AS  where 

A V„  =  -AS,  (20) 

kTlnlQ  v  1 

Equation  (20)  is  based  on  the  following  assumptions: 

1.  No  lateral  nonuniformities  (LNUs)  of  trapped  holes  in  the  oxide. 

2.  The  net  charge  of  interface  states  at  midgap  is  zero  and  are  uniform  between 
midgap  and  inversion.4 

3.  a  >  C/)/C0iC,  and  a  {Cd  +  C/r)  /Cor. 

The  assumption  that  a  »  {Cd  +  C/r)  /C0I  is  a  fairly  stringent  restraint.  As 
the  number  of  interface  states  increase,  C/r  also  increases.  This  assumption  essentially 
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limits  the  application  of  this  method.  At  high  total  dose  levels,  this  assumption  should 
be  checked. 

Similar  to  the  midgap  voltage  method,  the  voltage  shift  due  to  oxide  trapped 
holes,  A Vot>  is  calculated  using  the  shift  in  threshold  or  inversion  voltage  from  the 
relation 

A Vqt  =  A VTh  (AV/tfv-)  -  A VIT.  (21) 


2.3  AN w 


The  separation  of  AVmv  or  A Vth  into  AVot  and  AVjt  is  a  very  useful 
engineering  tool  for  the  evaluation  of  radiation  hardening  techniques  to  determine  the 
relative  amounts  of  hole  trapping  and  interface  state  generation.  However,  although 
a  reasonable  estimate  of  ANot  can  be  obtained  from  equation  (10),  a  reasonable 
estimate  of  A  Nit  is  not  necessarily  obtained  from  a  similar  relation  for  A  Vjr-  The 
value  of  ANit  should  be  determined  using  the  subthreshold  swing,  S,  as  described  by 
equation  (4). 

If  there  are  a  significant  number  of  interface  states,  the  capacitance  asso¬ 
ciated  with  these  states,  C/T,  is  in  parallel  with  Cp.  This  effects  the  value  of  the 
sub  threshold  swing,  S,  and  equation  (4)  becomes 


i + Cp+ c,t]  /ft  _4  (Cb+ 

Cqz  J  /  Or  \  C0z  / 


(22) 


Assuming  that  the  pre-irradiation  value  of  interface  state  density  is  small  (e.g.,  <  1010 
cm-2),  the  irradiation  induced  change  in  subthreshold  swing  is  given  by 
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where  Si  is  determined  using  equation  (4).  From  the  ratio  of  the  post- irradiation  to 
pre-irradiation  swing,  a  value  of  C/r  can  be  calculated  in  terms  of  5^,  Sj ,  Cp,  a  and 
Cox  as  follows. 


Cit  = 


(B~\ 

[tf  +  4|l 

(1  +  fc) 

1  - 
1  a2 

\Sx 
[  St 

(!  +  fe) 

1  -H 

}1/2> 

V 

(i  +  fc) 

_2_ 
a 2 

; 

cM  -  c 


D, 


(24) 


where 


H  =  1 


2 


The  depletion  capacitance  can  be  determined  from  the  pre-irradiation  swing  as  follows 


The  depletion  capacitance  defined  by  equation  (25)  includes  the  effects  of 
process  induced  interface  states.  Using  the  pre-  and  post-irradiation  subthreshold 
curves,  a  plot  of  the  energy  density  of  interface  states,  D/T,  versus  <f>s  can  be  deter¬ 
mined  in  a  manner  similar  to  the  Terman  technique16  for  MOS  capacitor  C-V  curves. 
If  the  values  of  Sn  and  5/  are  determined  at  several  values  of  ID  between  Img  and 
Iinvi  the  value  of  Dit  at  each  Ip  can  be  calculated  using  equations  (24)  and  (25), 
an4  the  relation  D!T  =  CiT/q.  The  total  number  of  interface  states  between  midgap 
and  inversion  is  found  by  integrating  Dit  over  energy  as  follows. 

r24>B 

ANit  =  /  DiTd(j>.  (26) 

If  Dit  is  constant  between  <f>s  =  <f>e  and  2<f>p ,  then 

A  Nit  =  (27) 

9 

Assuming,  a)  both  the  pre-irradiation  and  post-irradiation  subthreshold  channel  cur¬ 
rent  are  exponentially  dependent  on  Vp,  b)  the  irradiation  induced  interface  state 
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density 

CD/C0X 


is  constant  between  midgap  and  inversion,  and  c)  “a”  is  large  compared  to 
,  then  A  Nit  can  be  found  directly  from  AS  or  AVjt- 


A  Nit  = 


C0XAVit  =  A  SC0X<j>B 
q  kTtn\Q  ’ 


(28) 


where  AVjt  is  determined  from  equation  (14)  using  AViNV- 
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SECTION  3 


DISCUSSION  OF  VARIABLES  AND  ASSUMPTIONS  IN  THE 
SUBTHRESHOLD  CHARGE  SEPARATION  TECHNIQUE 


Analysis  of  a  MOSFET  that  has  a  near  linear  subthreshold  swing  using  the 
subthreshold  charge  separation  technique  is  well  defined.  The  subthreshold  swing  is 
extrapolated  to  find  the  midgap  voltage.  The  extrapolation  is  necessary  because  the 
midgap  current  is  below  the  leakage  of  the  device  and/or  the  instrumentation  is  not 
capable  of  measuring  such  low  current  levels.  When  a  device  exhibits  a  nonlinear 
subthreshold  swing,  one  must  choose  what  part  of  the  curve  to  extrapolate.  In  some 
cases,  the  choice  can  yield  widely  different  midgap  voltages,  which  affects  the  trapped 
hole  and  interface  state  calculations.  The  correct  choice  depends  on  the  purpose  of 
the  analysis  and  on  the  robustness  of  the  subthreshold  charge  separation  technique. 
There  are  inherent  assumptions  in  the  theoretical  development,  and  violations  of  these 
assumptions  may  or  may  not  greatly  affect  the  analysis  results.  To  analyze  nonlin¬ 
ear  characteristics,  the  analyst  must  hypothesize  possible  causes  and  determine  how 
the  subthreshold  technique  is  affected.  Combining  these  hypotheses  with  the  pur¬ 
pose  of  the  experimentation,  the  analyst  must  determine  the  portion  of  the  curve  to 
extrapolate. 


3.1  MOBILITY. 

Mobility  is  an  important  parameter  in  equation  (1).  Two  mobility  regions 
can  be  defined  in  a  MOSFET: 

1.  The  mobility  which  occurs  in  strong  inversion  (measured  via  the  transconduc¬ 
tance);  and 

2.  The  mobility  in  weak  inversion  or  subthreshold. 
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The  strong  inversion  mobility  is  greatly  affected  by  the  surface  of  the  silicon 
silicon-dioxide  interface,  since  the  inversion  channel  is  often  less  than  300  A.  In 
the  inversion  region,  the  vertical  electric  field  is  much  greater  than  the  transverse 
electric  field.  This  pulls  the  carriers  to  the  interface,  lowering  mobility  due  to  the 
roughness  (deviation  from  the  ideal  plane)  of  the  silicon/silicon  dioxide  interface.  The 
rougher  the  interface  the  lower  the  mobility.  In  the  subthreshold,  the  electric  field 
is  smaller,  and  the  interface  roughness  has  less  effect  on  the  mobility.  Thus,  as  the 
surface  potential  is  lowered  from  strong  to  weak  inversion,  the  interface  roughness 
has  less  effect  and  the  mobility  increases.  The  mobility  in  equation  (1)  is  the  weak 
inversion  mobility  and  can  be  substantially  different  than  the  strong  inversion  mobility 
determined  from  transconductance. 

Bulk  mobility  in  silicon  is  a  function  of  total  impurity  concentration  (the 
number  of  charged  impurity  atoms)  and  not  net  impurity  concentration,  which  deter¬ 
mines  the  carrier  concentrations.  When  both  N-  and  P-type  impurities  exist  in  silicon, 
acceptors  receive  electrons  from  the  donors  and  all  impurity  atoms  are  charged  (as¬ 
suming  total  ionization).  Thus,  in  theory,  the  weak  inversion  mobility  can  be  no  larger 
than  the  bulk  mobility. 

The  value  used  for  mobility  affects  the  results  of  the  subthreshold  charge 
separation  technique.  The  value  of  mobility  determines  midgap  and  inversion  currents 
and  hence  the  midgap  and  inversion  voltages.  Since  the  log  of  the  channel  current 
is  linear  with  gate  voltage,  errors  in  the  mobility  appear  as  the  log  of  the  error  on 
midgap  and  inversion  voltages,  unless  irregularities  exist  in  the  subthreshold  swing. 
Also,  since  the  trapped  hole  and  interface  state  voltage  shifts  are  calculated  from  the 
difference  between  the  pre-  and  post-irradiation  midgap  and  inversion  voltages,  the 
effect  of  an  error  in  mobility  is  further  decreased.  For  MOSFETs  with  near-linear 
subthreshold  characteristics,  the  technique  is  robust  to  the  actual  value  used.  If  a 
measured  threshold  voltage  is  used  instead  of  the  theoretical  inversion  voltage,  the 
error  from  an  incorrect  mobility  would  be  larger,  since  only  the  midgap  voltage  is 
affected  by  the  choice  of  the  mobility.  Irregularities  in  the  subthreshold  swing  could 
enhance  these  errors. 
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The  charge  separation  subthreshold  technique  assumes  that  the  mobility  is 
constant  from  ir'dgap  to  inversion, 


dlnfi 

d/3<f>s 


=  0. 


(29) 


If  this  derivative  is  nonzero,  the  subthreshold  swing  would  be  affected.  For  a  nonzero 
mobility  derivative,  equation  (4)  would  become 


dlnfi 

d(3<j>s 


(30) 


If  the  mobility  increases  as  surface  potential  is  lowered,  the  subthreshold 
swing  would  increase;  and  if  mobility  decreases,  the  subthreshold  swing  would  de¬ 
crease.  Fast  changes  in  mobility  over  short  ranges  in  surface  potential  can  cause 
distortions  in  the  subthreshold  swing.  For  example,  if  the  mobility  increases  from  a 
fixed  value  to  another  over  a  given  surface  potential,  then  the  subthreshold  swing  will 
increase  over  that  range  of  surface  potential.  The  subthreshold  swing  will  be  near 
identical  outside  that  range  as  shown  in  Figure  3. 

Channel  mobility  as  a  function  of  surface  potential  has  been  studied.  Mea¬ 
surements  of  mobility  have  not  been  done  for  surface  potentials  much  below  inversion. 
The  effect  of  surface  potential  and  radiation  induced  trapped  holes  and  interface  states 
on  channel  mobility  has  been  suggested.17  Figure  4  shows  a  possible  scenario  for  a  typi¬ 
cal  unirradiated  device.  Little  change  is  expected  to  occur  in  the  subthreshold  region. 
It  has  been  generally  accepted  that  ionizing  radiation-induced  interface  states  and 
trapped  charges  (holes)  within  50  A  of  the  surface  degrade  mobility.11  Interface  states 
act  as  charged  scattering  centers  and  degrade  mobility  similar  to  impurity  doping. 
Nonuniformity  of  interface  states  and/or  trapped  holes  can  create  further  degradation 
in  the  mobility  by  causing  minority-carrier  fluctuations,  which  have  a  large  effect  near 
inversion.  For  few  trapped  charges,  the  mobility  is  expected  to  decrease  with  surface 
potential  above  inversion,  but  for  many  trapped  charges,  it  can  increase  with  surface 
potential  because  of  the  decreasing  effect  of  minority  carrier  fluctuations.  Figure  5 
shows  possible  effects  of  the  radiation  induced  trapped  charges.  The  curves  of  inter¬ 
est  are  the  ones  which  show  mobility  increasing  as  the  surface  potential  approaches 
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Figure  3.  Effect  of  varying  mobility  on  MOSFET  subthreshold  curve. 


Figure  5.  Possible  effect  of  trapped  charges  on  mobility  for  an  irradiated 
MOSFET. 
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midgap.  The  hypothesis  for  this  occurrence  relies  on  the  assumption  that  most  inter¬ 
face  states  have  neutral  charge  at  midgap.  This  implies  the  post-irradiation  mobility 
near  midgap  should  improve  toward  the  pre-irradiation  condition.  Since  studies  have 
shown  post-irradiation  mobility  is  degraded  in  the  subthreshold,17  the  subsequent  in¬ 
crease  seems  possible,  due  to  the  decrease  in  charged  interface  states.  The  question 
is  how  rapidly  does  the  mobility  increase  as  the  surface  potential  approaches  midgap. 
The  increase  in  mobility  near  inversion  is  caused  by  a  decreased  effect  from  minority 
fluctuations. 

If  the  mobility  follows  curve  3  in  Figure  5,  the  effect  on  the  subthreshold 
swing  would  be  noticeable.  At  some  point  on  the  subthreshold  curve,  as  the  surface 
potential  approached  midgap,  the  swing  would  increase.  This  could  cause  a  second 
swing  or  increasing  swing  condition  as  is  shown  in  Figure  6.  The  lower  portion  of  the 
curve  has  a  larger  voltage  swing  per  decade.  In  this  case,  it  would  be  caused  by  a 
changing  mobility  and  not  additional  interface  states  in  the  bandgap.  Extrapolation  of 
measured  data  to  find  the  midgap  voltage  requires  judgement.  The  subthreshold  curve 
will  not  continue  at  the  larger  swing,  when  the  mobility  is  near  the  maximum  value. 
Extrapolating  from  the  top  or  bottom  of  the  curve  may  not  yield  the  correct  midgap 
voltage.  An  example  is  given  in  Figure  7  for  the  effect  of  the  increasing  mobility 
as  shown  in  Figure  6  on  the  subthreshold  curve.  To  analyze  this  data  using  the 
subthreshold  technique,  one  must  decide  from  which  part  of  the  curve  to  do  the  linear 
extrapolation.  Since  the  subthreshold  technique  assumes  that  mobility  is  constant  and 
uses  equation  (4),  the  first  swing  (top  part  of  the  curve  in  Figure  7)  should  be  used. 
This  will  give  the  voltage  shifts  caused  by  the  interface  states  and  trapped  holes  in 
the  oxide.  If  the  second  swing  is  used,  the  technique  will  give  the  voltage  shifts  caused 
by  the  number  of  interface  states  and  trapped  holes  necessary  to  produce  the  dashed 
line  from  midgap  to  inversion,  in  Figure  8.  Effectively,  this  procedure  is  attempting 
to  equate  a  number  of  interface  states  and  trapped  holes  to  produce  the  combined 
trapped  charge  and  mobility  effect  on  the  subthreshold  swing. 

Assuming  a  nonconstant  mobility,  the  choice  of  which  part  of  the  curve  to 
use  in  Figure  7  to  extrapolate  depends  on  the  purpose  of  the  analysis.  To  find  the 
correct  number  of  trapped  holes  and  interface  states  in  the  gate  oxide,  the  top  part 
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subthreshold  curve. 
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Figure  8.  Voltage  shifts  determined  from  second  slope  are  identical  to  the 
dashed  line  from  midgap  to  inversion. 


of  the  curve  should  be  used.  Extrapolating  the  lower  part  of  the  curve  would  give  a 
conservative  engineering  approximation  for  the  radiation  hardness  of  the  device,  based 
on  AVor  and  AVjt,  i.e.,  large  numbers  for  both  AVor  and  A V/y.  The  engineering 
approximated  total  dose  degraded  I-V  characteristic  of  the  device  would  be  equivalent 
to  the  dashed  line  from  midgap  to  inversion  in  Figure  8. 


3.2  CHANNEL  DOPING. 

Present  processing  methods  usually  will  yield  nonuniform  doping  in  MOS- 
FET  channel  regions.  Vertical  nonuniform  doping  occurs  in  MOSFETs  due  to  the 
threshold  adjust  implant.  This  implant  can  increase  the  subthreshold  swing  depend¬ 
ing  on  the  depth  of  the  implant,  but  it  will  not  cause  irregularities  in  the  swing. 
The  initial  increase  in  the  swing  does  not  affect  the  subthreshold  charge  separation 
technique  because  of  the  comparison  of  pre-  to  post-irradiation  curves. 

MOSFETs  which  have  nonuniform  doping  in  the  lateral  direction  under  the 
gate  cause  a  bending  of  the  conduction  and  valence  bands  with  respect  to  the  fermi 
level.  As  a  result,  the  entire  silicon  channel  will  not  be  at  the  same  potential.  When 
the  source  end  of  the  channel  is  at  midgap,  the  remainder  of  the  channel  will  be  above 
midgap.  It  is  necessary  to  determine  an  effective  channel  doping  to  use  in  equation 
(l).  Assuming  the  device  does  not  turn  on  until  the  entire  channel  is  inverted,  the 
effective  channel  doping  would  be  the  larger  value  near  the  source. 


3.3  LATERAL  NONUNIFORMITIES  (LNUs). 

In  recent  years,  LNUs  have  been  identified  as  a  cause  for  anomalous  results 
from  the  application  of  the  subthreshold  charge  separation  technique.10’1*  LNUs  ap¬ 
pear  to  cause  a  stretchout  of  the  subthreshold  swing  similar  to  the  stretchout  caused 
by  interface  states.  Consequently,  the  voltage  shift  due  to  interface  states  would  be 
overestimated  and  would  result  in  an  incorrect  radiation  characterization. 
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LNUs  refer  to  lateral  nonuniformities  of  trapped  holes  in  the  gate  oxide. 
The  surface  potential  at  the  silicon/silicon  dioxide  interface  is  directly  effected  by 
the  LNUs.  Irregularities  of  surface  potential  across  the  MOSFET  channel  can  cause 
preferred  current  paths,  so  at  a  given  surface  potential  the  channel  is  not  conducting 
current  uniformly.  The  effect  of  LNUs  on  a  MOSFET  channel  might  be  simulated  as 
many  parallel  transistors  all  with  different  threshold  voltages  but  identical  subthresh¬ 
old  swings.19  When  the  individual  curves  are  combined,  the  resulting  subthreshold 
swing  is  larger.  The  greater  the  spread  in  the  individual  threshold  voltages,  the 
greater  the  increase  in  the  subthreshold  swing. 

An  LNU  caused  increase  in  the  subthreshold  swing  is  interpreted  as  an 
increase  in  the  number  of  interface  states  by  the  subthreshold  charge  separation  tech¬ 
nique.  This  false  indication  of  interface  states  is  the  most  serious  effect  of  LNUs  and 
can  render  the  subthreshold  charge  separation  technique  as  an  inappropriate  analysis 
tool.  Lateral  nonuniformities  of  trapped  holes  violates  an  assumption  of  the  sub¬ 
threshold  charge  separation  technique  that  the  trapped  holes  are  uniform  throughout 
the  oxide  over  the  channel. 

It  appears19  that  the  creation  of  LNUs  are  radiation  bias,  dose  rate,  and 
energy  dependent.  Low  energy  high  dose  rate  radiation  pulses  such  as  from  an  x- 
ray  source  seem  to  produce  greater  LNU  effects.  Also,  the  effects  of  LNUs  appear 
to  anneal  at  room  temperature  in  relatively  short  times  (less  than  an  hour),19  but 
are  more  stable  at  cryogenic  temperatures.18  Annealing  of  LNUs  would  reduce  the 
subthreshold  swing  and  give  the  appearance  of  interface  state  annealing. 

LNUs  in  an  irradiated  MOSFET  could  possibly  be  identified  by  apparent 
annealing  of  interface  states  (reduction  of  subthreshold  swing)  at  room  temperature 
at  early  times  after  irradiation.  The  actual  number  of  interface  states  in  the  presence 
of  LNUs  can  be  estimated  using  the  subthreshold  technique.  If  LNUs  anneal  at  room 
temperature,  one  can  assume  that  the  LNUs  have  annealed  after  the  subthreshold 
swing  has  stabilized.  The  value  of  the  number  of  interface  states  calculated  using 
the  stable  subthreshold  swing  could  be  used  to  bound  the  true  value  at  earlier  times. 
Assuming  interface  states  are  relatively  stable  at  room  temperature,  this  value  of 


31 


the  number  of  interface  states  would  represent  a  worst  case  value  for  all  times  after 
irradiation  prior  to  this  measurement. 


3.4  INTERFACE  AND  TRAPPED  CHARGE  ASSUMPTIONS. 

Separation  of  radiation  induced  threshold  voltage  shifts  into  AVor  and  A  Vjy 
using  the  subthreshold  charge  separation  technique  is  sensitive  to  the  assumptions 
about  the  oxide  trapped  charge  and  interface  states.  The  subthreshold  technique 
assumes  that  the  interface  states  are  amphoteric.4  Only  acceptor  states  are  above 
midgap  and  donors  are  below  midgap.  This  provides  the  basis  for  assuming  that  all 
interface  traps  are  neutral  at  the  midgap  potential,  Vmg-  Assuming  that  only  holes 
are  trapped  in  the  oxide,  the  midgap  voltage  shift,  AVmG,  is  representative  of  the 
quantity  of  trapped  holes  in  the  oxide. 

Research  has  been  performed  which  proposes  the  trapping  of  electrons  in 
the  oxide20  or  at  the  interface21  in  hardened  oxides.  Neither  of  these  situations  would 
cause  non-linear  subthreshold  characteristics  to  occur.  The  trapped  electrons  would 
simply  have  a  cancelling  effect  on  the  oxide  trapped  holes.  In  this  case,  the  midgap 
voltage  shift  would  represent  the  difference  between  the  trapped  holes  and  electrons 
and  still  provide  a  good  engineering  tool  for  judging  radiation  hardness. 

Violation  of  the  assumption  that  interface  states  are  amphoteric  would  cause 
a  non-neutral  interface  charge  at  midgap,  but  would  not  affect  the  linearity  of  the 
subthreshold  characteristics.  Research  has  been  performed  showing  interface  states 
are  amphoteric4,22  and  nonamphoteric.23  In  n-channel  devices,  donors  above  midgap 
and  in  p-channel  devices,  acceptors  below  midgap  will  result  in  an  overestimation  by 
the  subthreshold  technique  of  the  number  of  trapped  holes  and  an  underestimation  of 
the  number  of  interface  states  by  the  number  of  nonamphoteric  interface  states. 
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3.5  PARASITIC  EFFECTS. 


Nonideal  subthreshold  characteristics  can  be  caused  by  parasitic  effects. 
These  effects  will  cause  aberrations  from  the  theoretical  linear  subthreshold.  Par¬ 
asitic  effects  can  include: 

1.  Junction  leakage,  drain  to  substrate,  and  gate  leakage  to  the  source  and/or 
drain. 

2.  Edge  and  parallel  parasitic  transistors. 

3.  Backchannel  and  edge  leakage,  and  kink  effect  in  a  SOS  or  SOI  MOSFET. 


3.5.1  Junction  Leakage. 

Leakage  currents  mask  the  actual  channel  subthreshold  current  below  the 
leakage  current  level  and  are  similar  to  instrumentation  measurement  noise  levels  in 
the  application  of  the  subthreshold  charge  separation  technique.  Only  the  part  of 
the  curve  which  is  above  the  leakage  is  used  for  the  extrapolation  to  midgap.  As 
with  leakage  current,  other  parasitic  effects  will  not  distort  the  subthreshold  until  the 
parasitic  current  becomes  near  equal  to  greater  than  the  MOSFET  channel  current. 
Parasitic  effects  are  not  related  to  the  hardness  of  the  gate  oxide  to  radiation-induced 
trapped  holes  and  interface  states,  which  the  subthreshold  charge  separation  technique 
calculates.  Parasitic  effects,  such  as  leakage,  may  be  related  to  the  radiation  hardness 
of  the  device  when  considered  in  conjunction  with  the  circuit  function. 

Leakage  currents  are  easily  distinguished  from  the  channel  subthreshold  cur¬ 
rent.  Subthreshold  analysis  of  this  data  has  always  been  done  using  the  data  above 
the  leakage.  Data  with  other  parasitic  effects  should  be  treated  similarly.  The  dif¬ 
ficulty  is  in  identifying  a  parasitic  effect.  Any  device  which  exhibits  a  subthreshold 
curve  with  significant  non-linear  distortions  could  have  a  parasitic  effect.  Applying 
the  subthreshold  charge  separation  technique  to  a  portion  of  the  curve,  which  contains 
a  parasitic  effect,  would  give  misleading  values  of  ANor  and  A  Nit-  Even  the  volt¬ 
age  shifts,  AFot  and  A V/r,  could  be  misleading,  since  these  values  imply  the  voltage 


33 


shifts  are  attributed  to  interface  states  and  trapped  holes  in  the  gate  oxide.  Again, 
this  reduces  to  using  the  subthreshold  charge  separation  technique  to  determine  the 
radiation  hardness  of  the  oxide,  or  the  midg^p  and  inversion  voltage  shifts  due  to  all 
the  effects  of  radiation  on  the  device. 


3.5.2  Edge  and  Parallel  Parasitic  Transistor  Effects. 

Edge  and  parallel  parasitic  transistor  effects  may  cause  distortion  in  the 
subthreshold  curve,  similar  to  what  has  been  observed  in  power  HEXFETs.  Figure  9 
shows  an  unirradiated  power  HEXFET,  which  has  a  bench  in  the  subthreshold.  This 
bench  has  been  observed  in  unirradiated  devices  from  different  manufacturers.  As  the 
total  dose  level  is  increased  for  these  devices,  the  pre-irradiation  bench  disappears, 
probably  due  to  the  dominance  of  the  radiation  effects.  The  power  HEXFET  may 
contain  10,000  or  more  individual  MOSFETs.  Nonuniformity  of  interface  channel 
doping  across  all  individual  devices  is  possible.  To  determine  the  effect  of  having 
a  few  MOSFETS  with  different  channel  interface  doping,  a  theoretical  subthreshold 
curve  was  generated  for  a  hypothetical  power  HEXFET.  The  model  was  run  for  100 
M  JSFETS  with  a  channel  interface  doping  of  1.0  X  1014  cm-3,  the  remaining  9,900 
had  an  interface  doping  of  1 .0  x  1016  cm-3.  This  model  gave  a  theoretical  subthreshold 
cjrve  shown  in  Figure  10.  The  upper  part  of  the  curve  is  identical  to  the  subthreshold 
swing  expected  for  the  bulk  of  the  devices.  The  lower  part  of  the  curve  is  dominated  by 
the  100  devices  with  lower  interface  doping.  Below  100  pA,  the  current  is  primarily 
through  the  100  devices  or  1.0  percent  of  the  entire  HEXFET.  The  bench  is  the 
transition  region  between  the  two  sets  of  MOSFETs.  The  drain  current  location 
and  voltage  width  of  the  bench  depends  on  the  percentage  of  the  nonuniformity  and 
difference  of  the  interface  channel  doping.  The  theoretical  model  has  a  bench  similar 
to  the  actual  power  HEXFET  in  Figure  9. 

A  subthreshold  bench  can  also  occur  in  a  discrete  MOSFET,  if  one  thinks  of 
a  single  MOSFET  as  many  MOSFETS  in  parallel.  If  an  edge  has  a  different  subthresh¬ 
old  curve  than  the  rest  of  the  device,  it  could  cause  a  bench  in  the  subthreshold.  Any 
nonuniform  dopin  5  defect,  which  forms  a  continuous  connection  between  the  source 
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Figure  9.  Power  HEXFET  pre-irradiation  subthreshold  curve. 


and  the  drain,  could  also  cause  a  bench.  As  was  shown  for  the  power  HEXFET,  these 
device  parasitics  can  affect  the  subthreshold  curve,  even  if  they  represent  only  1.0 
percent  of  the  device  width. 


3.5.3  Subthreshold  Effects  in  Insulated  Substrate  MOS  Transistors. 

The  subthreshold  characteristics  of  MOS  transistors  fabricated  on  insulating 
substrates  (i.e.,  silicon-on-sapphire  (SOS)  and  silicon-on-  insulator  (SOI))  are  affected 
in  three  ways.  These  include:  (1)  backside  leakage,  (2)  sidewall  leakage,  and  (3) 
floating  body  effects.  Backside  leakage  arises  from  a  parasitic,  a  conducting  path 
located  at  the  interface  between  the  epitaxial  silicon  and  the  insulator.  This  interface 
may  be  affected  by  charge  trapped  in  the  insulator  or  interface  states  existing  at  the 
boundary  between  the  two  materials.  If  the  source/drain  regions  extend  all  the  way 
to  the  silicon-insulator  boundary,  the  pre-irradiation  junction  leakage  is  increased  by 
the  recombination/generation  current  at  the  interface.  Post-irradiation  leakage  in  an 
n-channel  transistor  is  significantly  increased  because  charge  trapped  in  the  insulator 
tends  to  invert  the  silicon  near  the  boundary  and  creates  a  backside  channel  between 
the  drain  an'1  source.  The  result  is  that  the  lower  limit  of  subthreshold  leakage  current 
per  unit  width  is  initially  higher  in  insulated  substrate  devices.  In  post-irradiation  n- 
channel  devices,  the  drain- to-source  leakage  can  increase  to  the  point  that  it  obscures 
a  significant  portion  of  the  top-side  transistor  subthreshold  characteristic.  In  SOI 
technologies,  the  backside  leakage  path  can  be  turned  off  by  applying  a  negative 
potential  to  the  silicon  substrate. 

Sidewall  leakage  arises  from  the  parasitic  MOS  transistor  formed  at  the  ends 
of  the  intentional  MOSFET.  In  mesa  technologies,  the  gate  extends  over  the  edge 
of  the  silicon  island  and  forms  the  parasitic  transistor  on  a  different  silicon  plane. 
The  parasitic  device  may  have  a  lower  threshold  voltage  and  be  more  susceptible 
to  degradation  because  of  a  poorer  quality  gate  oxide.  A  similar  effect  is  found  in 
technologies  employing  a  trench  for  lateral  device  isolation.  A  parasitic  transistor 
is  formed  where  the  source  and  drain  abut  the  trench  sidewall.  The  trench  oxide 
is  subject  to  fringing  fields  from  the  drain-to-source  and  may  trap  charge  during 
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irradiation.  The  subthreshold  characteristic  of  both  mesa  and  trenched  devices  show 
changes  in  the  subthreshold  slope  as  a  function  of  gate  voltage.  The  turn-on  of  the 
parasitic  transistor  may  entirely  obscure  the  subthreshold  characteristic  of  the  top¬ 
side  transistor.  The  sidewall  transistor  is  eliminated  in  a  closed  geometry  transistor 
in  which  the  drain  completely  surrounds  the  gate  and  source  so  that  the  body  never 
encounters  the  sidewall. 

Floating  body  effects  arise  from  the  tendency  of  the  body  of  the  transistor 
to  develop  a  forward  bias  with  respect  to  the  source.  Leakage  or  avalanche  currents 
from  the  drain  to  the  body  give  rise  to  a  body  voltage  which  is  clamped  to  one  diode 
drop  above  the  source  voltage  by  the  source-to-body  diode.  The  result  is  to  move  the 
threshold  voltage  toward  depletion  through  the  body  effect  term  which  is  a  function  of 
the  square  root  of  the  body  voltage.  The  strong  inversion  I/V  characteristics  exhibit 
a  “kink”  effect  in  the  drain  current  versus  drain  voltage  curves.  The  subthreshold 
characteristics  are  influenced  by  the  shift  of  the  strong  inversion  threshold  voltage 
toward  depletion.  The  net  effect  is  to  produce  an  effective  change  in  slope  in  the 
portion  of  the  curve  near  high  injection  as  a  function  of  drain  voltage.  In  addition, 
subthreshold  curves  at  different  drain  voltages  tend  to  be  offset  toward  depletion  as  the 
magnitude  of  the  drain  voltage  increases.  The  floating  body  effect  can  be  mitigated 
by  providing  source-to-body  ties  as  part  of  the  process. 


3.6  CONCLUSIONS. 

The  overall  result  of  the  study  of  many  the  variables  which  affect  the  sub¬ 
threshold  curve  is  that  the  susceptibility  of  the  subthreshold  curve  to  these  variables 
increases  as  the  current  decreases  (excluding  the  kink  effect).  From  this  conclusion, 
extrapolation  of  nonideal  characteristics  (subthreshcld  characteristics  which  are  not 
described  by  equation  (l)),  which  occur  at  smaii  currents  (roughiy  below  100  pA),  is 
suspect.  The  nonideal  variation  may  represent  the  action  of  some  parasitic  or  defect 
in  the  device,  which  is  only  a  small  fraction  of  the  total  device.  If  this  is  true  and 
the  nonideal  variation  is  extrapolated,  the  subthreshold  charge  separation  technique 
would  give  results  which  are  not  representative  of  the  bulk  of  the  device.  The  midgap 


38 


voltage  determined  from  this  extrapolation  would  be  the  midgap  voltage  in  the  vicinity 
of  the  defect.  The  analysis  would  also  be  misleading  in  judging  the  radiation  hard¬ 
ness  of  the  gate  oxide.  From  an  engineering  view,  the  extrapolation  of  this  nonideal 
variation  still  would  not  be  justified.  Nonlinear  subthreshold  currents  at  or  below  100 
pA  would  not  usually  contribute  to  the  failure  of  a  circuit  function. 

To  minimize  the  effects  from  variables  which  cause  nonideal  subthreshold 
curves  at  low  currents,  the  upper  part  of  the  subthreshold  should  be  extrapolated. 
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SECTION  4 


APPLICATION 


In  order  to  apply  the  subthreshold  charge  separation  technique,  the  following 
data  set  and  parameters  must  be  known: 

1.  The  drain  and/or  source  current  (which  best  represents  the  channel  current) 
versus  gate  voltage  for  gate  voltages  corresponding  to  below  midgap  through 
strong  inversion.  This  data  set  may  be  taken  either  in  the  linear  or  saturated 
region  provided  VD  »  kT/q. 

2.  Measurement  of  mobility. 

3.  Gate  oxide  thickness. 

4.  Interface  channel  region  doping  density. 

4.1  MEASUREMENT  OF  ID(I0)  versus  VG. 

For  a  perfect  MOSFET,  the  drain  current  will  be  composed  almost  entirely 
of  the  channel  current  under  the  gate.  In  this  case,  the  drain  current  will  follow  an 
exponential  dependence  on  gate  voltage  in  the  subthreshold  with  a  single  valued  slope 
and  will  yield  ideal  MOSFET  subthreshold  characteristics.  In  real  MOSFETs,  there 
are  several  other  components  of  drain  cuirent  that  will  result  in  a  departure  from 
the  ideal  subthreshold  curve.  These  components  include  drain  to  substrate  leakage, 
gate  leakage,  edge  or  sidewall  leakage  and  back  channel  leakage.  The  edge  or  sidewall 
leakage  can  occur  under  a  field  oxide  or  an  extension  of  the  gate  oxide.  Both  sidewall 
and  back  channel  leakage  may  be  present  in  CMOS/SOS  or  CMOS/SOI  technolo¬ 
gies.  A  technique  for  reducing  the  interference  from  leakage  sources  is  to  measure  the 
source  current,  Is ,  rather  than  the  drain  current,  Ip.  This  will  eliminate  the  drain 
to  substrate  leakage  current,  since  source  and  substiate  are  usually  common,  and  it 
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may  reduce  other  leakage  components  as  well.  Before  a  decision  is  made  whether 
to  monitor  Ip  versus  VG,  or  Is  versus  VG,  the  currents  in  ail  available  leads  of  the 
MOSFET  should  be  monitored  as  a  function  of  VG  on  the  unirradiated  samples. 

The  test  configuration  for  the  measurement  of  channel  current  versus  gate 
voltage  is  shown  in  Figure  11  for  an  n-channel  MOSFET.  The  test  transistor  should 
have  all  four  regions  with  separate  contacts  if  possible.  For  many  technologies,  this 
may  not  be  possible.  For  example,  in  power  MOSFETs  the  source  and  substrate  are 
usually  shorted  with  metallization  on  the  chip.  In  CMOS/SOS  and  CMOS/SOI,  the 
substrate  is  often  floating.  In  these  cases  only  three  terminals  will  be  available.  If 
measurements  are  being  performed  on,  a)  test  transistors  located  on  process  control 
or  radiation  effects  test  chips,  or  b)  test  transistors  located  on  actual  microcircuit  die, 
the  test  transistors  should  be  designed  for  access  to  all  four  regions. 


4.1.1  Preliminary  Device  Characterization. 

In  order  to  determine  whether  the  drain  or  source  subthreshold  character¬ 
istics  best  represent  the  channel  current,  a  preliminary  test  should  be  conducted  on 
several  unirradiated  samples  to  determine  the  magnitude  of  the  leakage  components 
and  the  shape  of  the  subthreshold  curve.  If  all  four  regions  are  bonded  out  separately, 
then  the  drain  current,  Ip,  the  source  current,  Is ,  the  gate  current,  IG,  and  the  sub¬ 
strate  current,  IB,  should  be  measured  as  a  function  of  the  gate  voltage  from  below 
midgap  to  strong  inversion.  An  example  of  such  a  data  set  is  shown  in  Figure  12  for 
a  silicon  gate  n-channel  MOSFET  representative  of  a  3  /zm  bulk  CMOS  process.  In 
this  figure,  plotted  directly  from  the  graphic  output  of  the  HP4145,  all  four  terminal 
currents  are  shown  for  both  linear  (Vp  =  0.1  V)  and  saturated  ( Vp  =  10  V)  operation. 
The  top  trace  in  the  inversion  region  is  Ip  at  Vp  =  10  V.  Going  down  the  graph,  the 
second  curve  is  Is  at  Vp  =  10  V,  followed  by  Ip  at  VD  =  10  V,  The  fourth  curve  down 
is  both  Ip  and  Is  at  VD  =  0.1  V,  and  the  bottom  trace  is  composed  of  three  curves,  IB 
at  Vp  —  0.1  V  and  IG  at  both  Vp  =  10  V  and  0.1  V.  Because  of  the  large  vr  lue  of  Ip  at 
Vp  =  10  V,  the  source  current  is  a  much  better  representation  of  the  channel  current 
than  the  drain  current,  which  is  the  sum  of  channel  and  bulk  current.  Although  this 
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Figure  12.  Plot  of  all  four  terminal  currents  for  linear  and  saturated 
operation. 


difference  can  be  observed  in  pre-irradiation  curves  taken  in  saturation,  it  is  much 
more  pronounced  in  post-irradiation  curves  where  the  drain  to  substrate  leakage  has 
increased  due  to  the  surface  effects  around  the  periphery  of  the  drain.  If  bulk  or 
substrate  leakage  currents  are  appreciable,  the  source  current  is  a  better  measure  of 
the  channel  current.  However,  if  ‘here  is  any  appreciable  gate  leakage,  the  gate  may 
leak  to  the  source  causing  a  distortion  of  the  source  curve.  This  is  illustrated  in  Fig¬ 
ure  13,  which  shows  all  four  terminal  currents  at  Vd  =  0.1  V  and  5  V  for  a  p-channel 
device.  The  source  current  follows  the  gate  leakage  current  in  subthreshold,  whereas 
the  drain  current  is  a  good  measure  of  the  channel  current.  For  Vq  =  0  V  to  ~0.8  V, 
the  gate  current  and  bulk  current  are  equal  in  magnitude,  implying  that  in  this  range 
of  operation,  the  gate  leaks  to  the  bulk.  In  Figure  14,  the  four  terminal  currents  of  an 
n-channel  device  are  shown  illustrating  a  large  gate  leakage  current.  In  this  case,  the 
gate  leaks  to  the  substrate  and  the  source  and  drain  current  are  unaffected.  Thus,  if 
all  four  terminals  are  available,  I $  and  ip  are  near  ideal,  and  the  charge  separation 
technique  would  apply.  However,  if  this  were  a  three  terminal  device  with  the  source 
and  substrate  tied  together  internally,  I $  and  Id  would  no  longer  be  equal.  This  is  il¬ 
lustrated  in  Figure  15,  which  is  the  same  device  with  the  source  and  substrate  shorted 
together  at  the  device  pins.  In  this  case,  the  source  current  follows  the  gate  leakage 
current  until  the  channel  current  becomes  large  compared  to  the  gate  leakage.  For  a 
three  terminal  version  of  this  device,  the  drain  current  should  be  used  to  represent 
the  channel  current. 

The  above  examples  illustrate  how  the  gate  and  substrate  leakage  currents 
can  affect  the  source  and  drain  currents  in  subthreshold.  Therefore,  it  is  necessary  to 
fully  characterize  the  currents  in  all  available  device  leads  before  deciding  whether  to 
perform  data  analysis  on  the  source  current  or  drain  current. 


4.1.2  Pre-  and  Post-Irradiation  Measurement  of  Channel  Current 
Versus  Gate  Voltage. 

Once  the  preliminary  characterization  has  been  performed  to  determine 
whether  the  source  or  drain  current  is  more  representative  of  channel  current,  a  pre- 
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Figure  14.  Terminal  currents  of  an  n-channel  MOSFET  showing  gate 
leakage  to  substrate. 


******  GRAPHICS  PLOT  ****** 
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Figure  15.  Same  plot  as  previous  figure  with  substrate  and  source  shorted 
at  device  pins. 


irradiation  characteristic  is  measured.  The  measurement  of  Ip(Is )  versus  gate  voltage 
should  be  taken  from  below  midgap  to  strong  inversion.  The  drain  voltage  can  be  set 
for  either  linear  or  saturated  operation  in  inversion.  If  drain  voltage  is  set  for  lin¬ 
ear  operation,  the  condition  Vp  »  kT fq  must  be  satisfied.  At  room  temperature, 
Vp  should  be  at  least  100-150  mV.  For  saturated  operation,  Vp  can  be  set  at  the 
maximum  operating  voltage. 

There  are  several  precautions  that  should  be  considered  for  saturated  oper¬ 
ation: 

1.  Short  channel  effects  are  maximized  for  large  drain  voltage; 

2.  Hot  electron  effects  may  occur  for  large  drain  voltage; 

3.  Drain  to  substrate  currents  may  be  significant  for  large  drain  voltage.  This  effect 
is  illustrated  in  Figures  13  and  15; 

4.  In  floating  substrate  technologies  such  as  CMOS/SOS  and  CMOS/SOI  the  kink 
effect  may  occur  for  large  drain  voltage; 

5.  In  devices  with  Vt,:  >  Vp  [MAX)  saturated  operation  is  not  possible.  An  ex¬ 
ample  is  a  field  oxide  MOSFET  with  relatively  high  interface  doping  density; 
and 

6.  Measurement  of  mobility. 

Because  of  the  above  limitations,  it  is  recommended  that  the  subthreshold 
characteristic  be  measured  in  the  linear  region  with  100  mV<  Vp  <500  mV. 


4.2  MEASUREMENT  OF  MOBILITY. 

To  apply  the  subthreshold  charge  separation  technique,  a  value  for  the  chan¬ 
nel  mobility  is  required  and,  therefore,  a  measurement  of  the  transconductance  of 
the  device  is  needed.  The  transconductance  of  a  MOSFET  theoretically  depends  on 
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whether  the  device  is  in  the  linear  or  saturation  region.  In  the  linear  region,  the 
transconductance,  Gm,  is  given  by 


Gm  —  —nC0ZVD, 


and  in  the  saturation  region, 


Gm  —  —~nC0X  (VG  —  VTH  (Vim))  >  (32) 

where  m  is  a  function  of  doping  concentration  and  approaches  0.5  at  low  doping.  Both 
equations  (31)  and  (32)  could  be  used  to  determine  coefficients  for  the  drain  current 
of  equation  (1).  From  the  transconductance,  the  value  of  the  term,  K,  needed  for 
equation  (1)  is 

K  =  =  y^Goz  for  the  linear  region, 

Vo  L 


Gm  Z 

K  =  - — — - - — — - —  =  —fiC0Z  for  the  saturation  region.  (33) 

2m  [VG  —  Vth  (F;,vvj)  L 


For  the  linear  region,  one  just  divides  the  transconductance  by  the  applied 
drain  voltage  that  was  used  to  make  the  measurement. 

Transconductance  is  a  function  of  the  drain  and  gate  voltages  applied.  The 
drain  voltage  has  a  direct  effect,  but  also  has  an  indirect  effect  in  the  channel  length. 
The  mobility  is  a  strong  function  of  the  gate  voltage  due  to  the  vertical  electric  field,  as 
was  discussed  previously.  In  the  subthreshold  charge  separation  technique,  a  mobility 
is  needed  which  represents  the  weak  inversion  value.  Since,  for  a  MOSFET,  the 
mobility  increases  as  the  surface  potential  is  lowered,  the  maximum  strong  inversion 
value  should  be  used.  In  strong  inversion,  Figure  16  shows  the  mobility  decreasing 
as  gate  voltage  is  increased  for  a  3  nm  bulk  CMOS  technology  n-channel.  This  is 
an  important  consideration  in  making  linear  region  transconductance  measurements. 
Figure  17  shows  the  value,  K,  for  drain  voltages  of  0.1  and  5.0  volts.  When  the  drain 
is  at  5.0  volts,  the  value  of  K  is  decreased  by  almost  a  factor  of  3  which  is  caused  by 
the  decrease  in  mobility.  In  this  device,  the  effect  of  a  decreased  channel  length  is 
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Figure  16.  Plot  of  K  versus  gate  voltage  in  the  linear  region  for  drain 
voltages  of  0.1  and  1.0  volt. 


dominated  by  the  decrease  in  mobility.  This  implies  that  the  mobility  in  the  linear 
region  at  a  drain  voltage  of  5  volts  has  already  substantially  been  degraded  by  the 
gate  voltage.  Comparing  0.1  to  1.0  volt  in  Figure  16,  the  decrease  in  mobility  is  only 
about  15  percent.  For  this  process,  linear  region  mobility  degradation  due  to  the  gate 
electric  field  does  not  sharply  increase  until  the  drain  voltage  is  larger  than  1.0  volt. 


Saturation  transconductance,  equation  (32),  can  also  be  used  to  evaluate  the 
coefficient  for  the  drain  current.  This  technique  does  not  seem  to  yield  the  maximum 
values  that  are  obtained  in  the  linear  region.  Figure  18  shows  the  use  of  equation 
(32)  for  the  same  device  used  to  generate  Figures  16  and  17.  The  mobility  term  is  30 
percent  lower  when  determined  using  the  saturation  region. 

Mobility  is  needed  to  calculate  the  L.idgap  and  inversion  currents.  Since 
interface  states  degrade  the  mobility,  it  is  expected  that  post-irradiation  mobility  will 
be  less  than  the  pre-irradiation  value.  The  subthreshold  charge  separation  technique 
has  been  done  using  the  pre-irradiation  mobility  to  analyze  the  post-irradiation  data, 
and  recalculating  the  mobility  for  each  post-irradiation  data  set.  If  the  pre-irradiation 
value  is  used,  one  assumes  that:  1)  the  weak  inversion  mobility  changes  little  with 
irradiation,  or  2)  the  mobility  calculated  from  the  post-irradiation  curves  are  a  worse 
estimate  of  the  post-irradiation  weak  inversion  mobility  than  the  pre-irradiation  cal¬ 
culated  mobility. 

Using  the  interface  state  neutrality  assumption,  the  mobility  at  midgap  is 
better  approximated  by  the  pre-irradiation  value.  At  midgap,  all  interface  states 
located  above  midgap  have  neutral  charge,  and  it  is  assumed  that  there  are  not  many 
trapped  holes  near  the  interface  to  affect  the  mobility. 

Thus,  for  the  calculation  of  the  midgap  current,  the  nondegraded  pre-irradiation 
value  should  be  used,  for  post-irradiation,  the  mobility  is  assumed  to  be  degraded 
by  interface  states,  and  at  the  inversion  voltage,  it  is  also  assumed  that  all  interface 
states  are  charged.  Thus,  the  value  of  mobility  at  the  post-irradiation  inversion  voltage 
may  be  better  represented  by  the  mobility  calculated  from  that  data  set.  If  a  mea¬ 
sured  threshold  voltage  is  used,  only  the  pre-irradiation  transconductance  needs  to  be 
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measured.  If  the  inversion  current  is  calculated,  then  both  pre-  and  post-irradiation 
mobilities  could  be  measured  and  used  in  the  data  analysis. 

To  achieve  the  best  estimates  for  the  mobilities  at  midgap  and  post-irradiation 
inversion,  the  smallest  possible  drain  voltage  should  be  used.  This  enables  the  mobility 
to  be  measured  at  small  gate  voltages  and,  thereby,  minimizes  the  surface  roughness 
and  electric  field  effect.  Also,  the  mobility  would  be  measured  at  a  gate  voltage  close 
to  inversion  providing  a  better  estimate  of  the  actual  weak  inversion  mobility.  For 
pre-irradiation  data  where  interface  state  concentration  is  small,  the  mobility  should 
be  degraded  very  little  at  small  gate  voltage.  If  this  degradation  is  small,  the  measured 
mobility  should  be  close  to  the  actual  mobility  at  midgap. 


4.3  DATA  ANALYSIS. 

As  stated  earlier,  the  subthreshold  charge  separation  technique  can  be  used 
either  as  an  engineering  tool  to  compare  the  threshold  voltage  shift  contributions  from 
trapped  holes,  AVori  and  interface  states,  A V/r,  or  it  can  be  used  as  a  research  tool 
to  determine  the  energy  density  of  the  interface  states,  Djt .  For  use  as  an  engineering 
tool  to  either  characterize  technologies,  investigate  process  hardening  techniques  or 
perform  hardness  assurance,  the  data  analysis  is  quite  straightforward.  However,  if 
the  technique  is  to  be  used  as  a  research  tool  to  explore  the  spatial  dependence  of 
trapped  holes  or  the  energy  dependence  of  the  interface  states,  the  data  requirements 
are  more  severe  and  the  analysis  more  complicated. 


4.3.1  Separation  of  A Vth  into  AVOT  and  A V/r. 

The  data  analysis  procedure  for  the  separation  of  A  Vth  into  the  components 
AVqt  and  A V/r,  is  based  on  the  following  equations  as  discussed  in  Section  2, 


A  Vqt  =  AFjvrcj 
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and 


AV/r  =  AVy//  —  A  Vqti 


or 


and 


A  Vqt  —  A  Vth  —  A  Vjt- 


Using  the  midgap  method,  the  midgap  voltage  shift,  A Vmg,  is  calculated  from  the  pre- 
and  post-irradiation  subthreshold  log  Id  versus  Vq  characteristics  using  the  following 
procedure: 


1.  1)  Determine  the  value  of  f  fiC0X  =  ^  from  a  calculation  of  peak  transconduc¬ 
tance,  A/xj/AUg  using  the  pre-irradiation  Id  versus  Vg  data  in  strong  inversion 
as  shown  previously  in  Figure  16. 

2.  Obtain  values  of  the  oxide  thickness,  tox,  and  the  channel  (substrate)  doping, 
Na  or  Nd- 

3.  Calculate  the  channel  midgap  current,  Img,  from  equation  (12)  and  (l)  of  Section 

2. 

4.  Extrapolate  the  straight  line  pre- irradiation  subthreshold  log  Id  versus  Vq  char¬ 
acteristic  to  a  value  Img •  The  gate  voltage  at  Img  Is  Vmgi- 

5.  Extrapolate  the  straight  line  post-irradiation  subthreshold  log  Id  versus  Va  char¬ 
acteristic  to  the  same  value  Img-  This  gate  voltage  is  Vmgt 

6.  A Vmg  =  Vmg~i  ~  Vmgi  =  AVor- 

7.  Repeat  5)  and  6)  for  each  post-irradiation  subthreshold  characteristic. 


Because  the  oxide  trapped  holes  due  to  irradiation  is  positive,  AVer  will 
always  be  negative. 
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Using  the  subthreshold  swing  method,  the  voltage  shift  due  to  interface 
states,  AV/r  is  calculated  from  the  pre-  and  post-irradiation  log  Id  versus  VG  charac¬ 
teristics  using  the  following  procedure: 

1.  Obtain  the  value  of  channel  (substrate)  doping,  N \  or  ND. 

2.  Calculate  the  straight  line  fit  to  the  pre-irradiation  log  Id  versus  Vq  subthreshold 
characteristics. 

3.  Calculate  straight  line  fit  to  the  post-irradiation  Log  Id  versus  VG  subthreshold 
characteristics. 

4.  A S  =  S7  -  Si. 

5.  AV/r  =  ( q<f>B/kTtnlO)AS 

6.  Repeat  3),  4),  and  5)  for  each  post-irradiation  subthreshold  characteristic. 

In  order  to  calculate  the  value  of  AV/r  using  the  midgap  method,  or  AV0r 
using  the  subthreshold  swing  method,  the  irradiation  induced  threshold  voltage  shift 
A VTH  or  the  inversion  voltage  shift,  A Vmv  must  be  determined.  The  value  of  AVTH 
is  simply 

AVp//  =  VjH'i  ~  Vthi, 

where  Vthi  is  the  pre-irradiation  value  of  Vth  and  Vr//7  is  the  post-irT  idiation  value 
of  Vth •  The  threshold  voltage  is  the  extrapolated  gate  voltage  for  Id  —  0  using  either 
a  plot  of  \[Td  versus  Vq  for  operation  in  saturation  or  a  plot  of  Id  versus  Vq  for 
operation  in  the  linear  region.  The  calculation  of  Vth  is  illustrated  in  Figures  19  and 
20. 


The  value  of  AV!Nv  is  calculated  using  equation  (7)  and  the  following  pro¬ 
cedure: 

1.  Determine  the  value  of  ffiC0X  =  ^  from  a  calculation  of  peak  transconduc¬ 
tance,  AId/AVq  using  the  pre-irradiation  Id  versus  VG  data  in  strong  inversion 
as  shown  previously  in  Figure  16. 
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Threshold  voltage  determination  from  plot  of  y/lp  versus  VG 
for  saturated  operation  with  Vjj  —  5  V. 


2.  Obtain  values  of  the  oxide  thickness,  tox,  and  the  channel  (substrate)  doping, 
Na  or  ND. 

3.  Calculate  the  channel  inversion  current,  Iinv,  from  equation  (7). 

The  value  of  AVwv  is 

A Vjyy  =  Vjjvvt  ~  Vjhvi . 

where  Vwvi  is  the  pre-irradiation  value  of  Vjsv  and  V/jvv-y  is  the  post-irradiation 
value. 


In  p-channel  devices,  A Vth  will  always  be  negative  since  the  net  charge  of 
the  interface  states  is  positive,  causing  AV/r  to  be  negative.  In  n-channel  devices, 
where  the  net  charge  of  the  interface  states  is  negative,  AV/r  will  be  positive.  There¬ 
fore,  in  n-channel  devices,  A  Vth  may  be  positive  or  negative  depending  on  whether 
lAV/rl  is  larger  or  smaller  than  |AV0r|- 

If  the  midgap  method  is  used,  AV/r  is  calculated  by 
AV/r  =  i\.  rit  (AV/jvv)  -  AVor- 

If  the  subthreshold  swing  method  is  used,  AVor  is  calculated  by 

AV0r  =  A  Vth  (AV//yy)  —  AV/r. 

Figure  21  is  an  illustration  of  a  typical  pre-  and  post-irradiation  plot  of 
log  Ip  versus  Vb  for  an  n-channel  MOSFET  showing  the  separation  of  AVr/r  into 
AVor  ana  AV/r.  In  Figure  22,  a  p-channel  pre-  and  post-irradiation  subthreshold 
characteristic  are  shown. 

Although  the  midgap  method  of  the  subthreshold  voltage  separation  tech¬ 
nique  can  be  applied  graphically  as  shown  in  the  examples,  for  large  data  sets,  a 
computer  program  should  be  written  to  perform  the  necessary  calculations.  The  only 
input  data  required  for  such  a  program  are  the  (Id,  Vq)  point  pairs,  the  drain  voltage, 
Vd,  and  the  parameters  tox  and  N a  or  No  for  the  calculation  of  midgap  and  inversion 
voltages. 
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Figure  22.  Typical  voltage  shift  separation  for  p-channel  device. 


4.3.2  Use  of  Subthreshold  Data  for  Determination  of  A  Not  and  A  iV/y. 


The  subthreshold  technique  can  also  be  used  for  the  investigation  of  irradi¬ 
ation  induced  oxide  trapped  charge,  AJVor,  and  the  energy  density  of  the  interface 
states,  Dit,  or  its  integrated  value  from  midgap  to  inversion,  ANjt.  The  value  of 
A  Not  can  be  obtained  from  A  Vmg  provided  certain  assumptions  can  be  made  about 
the  spatial  distribution  of  the  trapped  holes.  A  value  of  Dir  versus  gate  voltage  (sur¬ 
face  potential)  can  be  calculated  provided  certain  assumptions  are  made  concerning 
the  subthreshold  mobility. 

If  the  net  charge  of  the  interface  states  at  midgap  is  zero,4  the  voltage  shift, 
AVmg,  is  due  entirely  to  the  irradiation  induced  trapped  hole  density  in  the  gate 
oxide.  Therefore, 

—  1  r^ox 

AVmg  =  —  /  A  p0X(x)xdx  (34) 

where  Ap0I(x)  is  the  spatial  positive  charge  density  in  the  oxide,  tox  the  oxide  thickness 
and  x  is  measured  from  the  silicon/silicon  dioxide  interface  to  the  gate.  If  A pox  is  a 
sheet  charge  next  to  the  silicon  containing  qANor  charges,  AVmG  =  qANor /Cox ;  and 
if  A pox  is  a  constant  equal  to  qANor,  AVmG  =  qANor  /2C0X. 

For  many  gate  oxide  MOSFETs  the  trapped  noles  after  hole  transport  and 
trapping  are  located  within  50-100  A  of  the  silicon-silicon  dioxide  interface.14  If  this 
condition  exists,  then  for  tox  »  50  -  100A, 

ANot  *  (35) 

9 

Unfortunately,  many  state-of-the-art  hardened  oxides  are  on  the  order  of  200-500  A 
thick.  In  this  case,  for  equation  (35)  to  be  valid  the  trapped  holes  would  need  to 
be  located  on  the  order  of  50  A  or  less  from  the  interface.  Since  hole  detrapping  by 
tunneling  begins  immediately  after  hole  transport  and  tiapping,  those  holes  trapped 
within  a  few  10s  of  A  of  the  interface  are  detrapped  in  times  comparable  to  the  post- 
irradiation  measurement  time.  Therefore,  in  gate  oxides  of  <^o00  A,  equation  (35)  may 
severely  underestimate  the  number  of  trapped  holes.  However,  since  the  assumption 


62 


of  uniform  trapping  yields  a  value  twice  the  value  for  Si-Si02  interface  trapping,  the 
underestimate  will  be  less  than  a  factor  of  2. 

Because  most  techniques  applied  to  MOS  capacitors  use  the  same  assump¬ 
tions  about  the  spatial  dependence  of  the  trapped  holes,  the  use  of  equation  (35)  will 
be  quite  adequate  for  comparing  ANor  to  values  obtained  from  other  techniques  for 
thin  gate  oxides. 

There  is  some  evidence15  that  in  thick  field  oxides,  including  bipolar  recessed 
field  oxides,  the  holes  are  uniformly  trapped  in  the  bulk.  In  this  case, 

AJVor  =  2AVMaC“'  (36) 

<7 

Equation  (36)  is  also  valid  for  irradiation  at  low  temperature  (<76°  K)  since, 
in  this  case,  hole  transport  is  impeded,  and  the  holes  are  uniformly  trapped  in  the 
bulk  near  where  they  are  created. 

The  calculation  of  A Njt  is  given  in  Section  2.3. 


4.4  APPLICATION  RECOMMENDATIONS. 

Use  of  the  midgap  or  subthreshold  swing  method  of  the  subtb  ,-shold  charge 
separation  technique  requires  a  straight  line  approximation  to  the  subthreshold  curve. 
The  midgap  method  uses  the  straight  line  fit  for  extrapolation  to  a  midgap  voltage, 
and  the  subthreshold  swing  method  uses  the  fit  to  calculate  the  voltage  shift  due 
to  interface  states.  The  midgap  method  has  been  more  widely  used  by  the  techni¬ 
cal  community.  Also,  the  additional  assumption  that  a  >  ( CD  +  C/j)  / Coz  hinders 
the  application  of  the  subthreshold  swing  method  to  technologies  which  exhibit  high 
interface  state  generation. 

Linear  fits  to  the  subthreshold  curve  is  complicated  by  nonideal  subthresh¬ 
old  characteristics.  The  approach  to  the  application  of  the  subthreshold  technique  to 
MOSFETs  with  nonideal  characteristics  is  discussed  in  Section  3.  Once  a  linear  fit  to 
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a  nonideal  subthreshold  curve  is  chosen,  the  analysis  follows  the  approach  presented  in 
Section  4.3.  The  calculated  voltage  shifts  from  MOSFETs  with  nonideal  characteris¬ 
tics  should  be  combined  with  speculation  to  the  cause  of  the  nonideal  characteristics. 
This  will  aid  other  readers  in  understanding  the  analysis. 
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APPENDIX 


RECOMMENDED  TEST  PROCEDURE  FOR  SEPARATING  AVOT 
AND  A V1T  FROM  AVTH  IN  IRRADIATED  MOSFETS  USING 

SUBTHRESHOLD 


MIL-STD-883 
TEST  METHOD  10XX 

Subthreshold  Technique  for  separation  of  ionization  induced  MOSFET  thresh 
old  voltage  shifts  into  voltage  components  resulting  from  oxide  trapped  charge  and 
interface  states. 


A.l  PURPOSE. 

The  purpose  of  this  test  method  is  to  separate  the  ionization  induced  thresh¬ 
old  voltage  shift,  A Vth,  in  an  MOS  transistor  into  components  resulting  from  oxide 
trapped  charge,  A Vqt,  and  irradiation  induced  interface  traps,  AVIT.  The  technique 
used  in  this  test  method  is  based  on  the  pre-  and  post-irradiation  electrical  character¬ 
istics  in  the  subthreshold  region. 


A. 1.1  Definitions. 

a.  Threshold  Voltage,  Vth •  The  value  of  gate  to  source  voltage  extrapolated 
to  a  drain  current,  Id,  of  zero  using  a  plot  of  Id  or  \JTd,  versus  gate  to  source 
voltage,  Vgs-  For  linear  region  operation,  Id,  is  extrapolated  and  for  saturated 
operation,  \JTd  is  extrapolated. 

b.  Radiation  Induced  Threshold  Voltage  Shift,  A  Vth-  The  post-irradiation 
value  of  Vth  minus  the  pre-irradiation  value.  A  Vth,  can  be  either  positive  or 
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negative  depending  on  whether  the  net  radiation  induced  charge  is  negative  or 
positive,  respectively. 

c.  Inversion  Voltage,  V/^rv  The  gate  to  source  voltage  corresponding  to  a  sur¬ 
face  potential  of  2 <f>B,  where  <f>s  is  the  bulk  Fermi  potential. 

d.  Radiation  Induced  Inversion  Voltage  Shift,  A V/atv*  The  post-irradiaticn 
value  of  Vim  minus  the  pre-irradiation  value. 

e.  Inversion  Current,  J/jvv.  The  value  of  drain  current  at  a  surface  potential  of 
2 

f.  Midgap  Voltage,  Vmg.  The  gate  to  source  voltage  for  a  surface  potential  of 
4>B‘ 

g.  Radiation  Induced  Midgap  Voltage  Shift,  A Vmq.  The  post-irradiation 
value  of  Vmg  minus  the  pre-irradiation  value. 

h.  Midgap  Current,  I  mg-  The  value  of  drain  current  at  a  surface  potential  of 
<I>B- 

i.  Surface  Potential,  <f>s .  The  potential  at  the  semiconductor  surface  measured 
with  respect  to  the  intrinsic  Fermi  level.  For  the  flat  band  condition,  <f>s  =  0. 

j.  Bulk  Fermi  Potential,  <f>B-  The  surface  potential  at  which  the  Fermi  level 
and  the  intrinsic  Fermi  level  coincide. 

k.  Transconductance,  The  rate  of  change  in  Id  with  respect  to  source  to 
gate  voltage  at  a  constant  source  to  drain  voltage.  For  an  ideal  MOSFET  in  the 
linear  region,  Gm  =  [Z/ Le/j)  hC0XVd,  where  Z  is  the  channel  width,  Le/f  the 
effective  channel  length,  n  the  mobility,  Cox  the  oxide  capacitance  per  unit  area 
and  V o  the  drain  voltage. 

l.  Oxide  Thickness,  toz.  The  thickness  of  the  gate  oxide  measured  from  the 
gate-Si02  interface  to  the  Si-Sify  interface. 

m.  Channel  Doping  Density,  Nd,Na.  The  donor  (P-type)  or  acceptor  (N-type) 
concentration  .  the  channel  region  adjacent  to  the  gate  oxide-silicon  interface. 

n.  Subthreshold.  <f>B  <  <f>s  <  2<j>B,  also  referred  to  *s  weak  inversion. 
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o.  Subthreshold  Swing,  S.  The  subthreshold  swing  is  the  reciprocal  of  the  sub- 
threshold  slope  of  log  Ip  versus  Vas^  S  is  measured  in  units  of  volts  per  decade. 

A. 2  APPARATUS. 

A. 2.1  Instrumentation. 

To  measure  the  subthreshold  drain  current  versus  the  gate  to  source  voltage 
on  standard  microcircuit  MOSFETS,  the  instrumentation  required  consists  of,  as  a 
minimum,  two  voltage  sources  and  four  ammeters  with  the  following  characteristics. 

Equipment  Requirements 

Drain-Source  Power  Supply  ±10  V,  0.01  V  resolution 
Gate-Source  Power  Supply  ±10  V,  0.001  V  resolution 
Ammeters  for  measuring  ±1  pA  minimum  resolution 
terminal  currents  ±10  mA  maximum  range 

If  the  MOSFETs  under  test  (MUTs)  are  power  MOS,  parasitic  field  oxide  structures 
-it  higher  voltage  linear  MOSFETs,  the  maximum  voltage  requirements  will  be  much 
higher.  For  example,  parasitic  field  oxide  transistors,  FOXFETs,  may  have  initial 
threshold  voltages  of  several  hundred  volts. 

In  most  applications  of  the  charge  separation  technique  the  subthreshold 
drain  current  versus  gate-source  voltage  characteristic  will  be  measured  with  a  pro¬ 
grammable  tester  with  the  proper  current  and  voltage  capabilities.  The  most  widely 
used  instrument  for  MOSFET  subthreshold  I-V  measurements  is  the  HP4145  param¬ 
eter  analyzer.  The  basic  system  has  a  100  V,  100  mA  capability  with  1  mV,  1  pA 
resolution.  For  rapid  subthreshold  I-V  measurements,  a  programmable  tester,  the 
Transient  Annealing  Test  System  (TATS)  has  been  developed  by  Design  Engineering, 
Inc.  (DEI),  of  Albuquerque,  NM.  This  instrument  is  current  programmable  and  is 
capable  of  100  pA  measurements  in  200  /is  and  >100  nA  measurements  in  70  /is. 
The  subthreshold  I-V  measurements  can  be  made  with  any  instrument  having  the 
appropriate  characteristics. 
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A. 2. 2  Temperature  Control. 


The  subthreshold  drain  current  varies  as  the  exponential  of  •  There¬ 
fore,  the  temperature  of  the  measurement  should  be  controlled  to  within  a  few  °C,  since 
the  chaise  extraction  requires  a  comparison  of  pre-irra  ”  i;on  and  post-irradiation 
data.  For  many  applications,  the  charge  separation  will  be  performed  at  a  temper¬ 
ature  other  than  room  ambient.  For  these  measurements,  special  heating  and/or 
cooling  apparatus  are  required.  The  heating/cooling  apparatus  snail  be  capable  of 
controlling  the  MUT  temperature  to  within  ±2°C.  A  temperature  sensing  device  shall 
be  provided  which  shall  be  in  thermal  contact  with  the  MUT.  The  readout  shall  be 
capable  of  0.5  C  resolution  and  calibrated  to  ±1°C  accuracy.  The  temperature  control 
apparatus  shall  allow  for  power  and  sense  leads  to  the  MUT  and  shall  not  degrade 
the  drain  current  measurement  by  more  than  10  pA. 


A. 3  PROCEDURE. 

A.3.1  Device  Identification 

in  all  '.-»ses,  each  MUT  shall  be  individually  identified  (marked  with  a  per¬ 
manent  marking)  and  its  identification  shall  be  included  in  all  pre-  and  post-irradiation 
data. 


A. 3. 2  Initial  Characterization. 

Using  the  appropriate  test  apparatus,  the  MUT  shall  be  characterized  using 
the  test  configuration  illustrated  in  Figure  23.  The  gate-source  voltage  shall  be  swept 
from  accumulation  to  inversion  in  order  to  clearly  identify  the  leakage,  subthreshold, 
and  inversion  regions.  The  drain  voltage  may  be  selected  for  either  linear  or  satu¬ 
rated  operation.  It  is  recommended  that  measurements  be  made  in  the  linear  region 
with  Vd  between  100  mV  and  150  mV.  As  the  gate-source  voltage  is  stepped  from 
leakage  (accumulation)  tJr'n.gh  inver;  on,  all  four  currents,  1g,Is,Id  and  shall 
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be  monitored.  If  there  is  no  access  to  the  body  (channel,  substrate)  then  only  three 
currents  shall  be  monitored.  In  some  devices  (SOS,  SOI),  the  body  floats  and  in  some 
devices  (power  MOS)  the  source  and  body  are  internally  tied  together.  Once  all  of 
the  available  terminal  currents  are  measured,  a  plot  is  made  of  log  I  versus  VGs-  All 
three  or  four  currents  should  be  plotted  on  the  same  graph  for  comparison. 

The  purpose  of  plotting  all  currents  is  to  analyze  the  MUT  for  leakage  paths 
and  parasitics  to  determine  the  idealty  of  the  subthreshold  characteristic.  Discussion 
of  nonideal  behavior  is  given  in  the  guideline  document.  Distortions  of  the  subthresh¬ 
old  characteristic  may  be  caused  by  gate  leakage,  body-drain  or  body-source  leakage, 
or  parasitics.  After  pre-irradiation  characterization,  it  may  be  decided  that  the  part 
is  unacceptable  for  charge  separation  studies  or  that  the  source  current  rather  than 
drain  current  is  a  better  representation  of  the  ideal  channel  current. 

An  example  of  a  near  ideal  subthreshold  characteristic  is  shown  in  Figure  24 
for  an  n-cham  el  transistor.  For  this  example,  the  leakages  are  all  small  and  the  source 
or  drain  current  could  be  used  for  charge  separation. 

Based  on  the  initial  characterization,  a  decision  must  be  made  as  to  the 
suitability  of  the  test  sample  for  charge  separation.  For  those  samples  deemed  suitable, 
a  decision  is  made  whether  to  measure  drain  current,  Id  or  source  current,  Is-  If  both 
Id  and  Is  are  suitable,  then  the  decision  may  be  based  on  convenience. 


A.4  PRE- /POST-IRRADIATION  CHARACTERIZATION. 

A.4.1  Pre-Irradiation  Characterization. 

Because  of  the  sensitivity  of  Id(s)  temperature,  the  pre-irradiation  charac¬ 
terization  used  for  charge  separation  should  be  made  either  in  a  temperature  controlled 
chamber  or  immediately  prior  to  irradiation.  The  drain  (or  source)  current  is  mea¬ 
sured  from  accumulation  (or  leakage)  through  the  subthreshold  region  and  well  into 
inversion  at  a  fixed  value  of  drain  to  source  voltage  | [Vds ).  For  linear  measurements, 
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decade/div 


ID, (A)  IS, (A) 


Figure  24.  Near  ideal  subthreshold  characteristics  from  an  n-channel 
transistor. 
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decade/div 


VDS  should  be  chosen  between  100  and  150  mV.  For  saturated  operation,  the  value 
of  VDs  will  depend  on  the  threshold  voltage  and  the  rated  value  of  Vds ■  If  the  test 
samples  are  conventional  CMOS  microcircuit  MOSFETs  and  saturated  measurements 
are  specified,  then  VDs  can  be  set  at  the  nominal  power  supply  voltage. 

If  the  measuring  instrument  is  automatically  stepped  through  programmed 
values  of  VG5,  the  step  size  shall  be  adjusted  so  that  a  minimum  of  3  Id{s)  values  per 
decade  current  are  measured  in  subthreshold.  Also,  the  dwell  time  at  a  fixed  Vcs  shall 
be  adjusted  to  assure  adequate  settling  time.  If  an  HP4145  or  similar  instrument  is 
used,  the  medium  or  long  measurement  time  setting  should  be  used. 

The  value  of  Id(s)  for  each  value  of  shall  be  recorded,  along  with  the  part 
identification  number,  date  and  time,  temperature  and  a  notation  for  pre-irradiation. 


A.4.1.1  In-Source  and  In-Situ  Test  Setup.  If  the  MUT  is  being  measured  in¬ 
source  (i.e.,  while  it  is  being  irradiated),  or  in-situ  (i.e.,  in  the  irradiation  test  fixture 
but  not  while  it  is  being  irradiated),  the  test  fixture  and  cabling  connecting  the  test 
instrument  must  be  designed  for  low  current  measurements.  In  addition,  the  test 
fixture  must  provide  a  mechanism  for  switching  the  MUT  from  the  irradiation  and/or 
annealing  bias  conditions  to  the  measurement  mode.  In  order  to  obtain  accurate 
measurements  in  the  subthreshold  region,  the  leakage  and  noise  currents  in  the  test 
fixture  and  cabling  must  be  less  than  10  pA.  Therefore,  measurement  cables  should 
be  shielded  and  of  minimum  length.  The  leakage  and  noise  currents  shall  be  measured 
in  the  test  setup  by  monitoring  the  Ip(S)  versus  VGS  characteristic  v  'th  no  device  in 
the  socket.  Maximum  current  over  the  voltage  range  of  interest  should  not  exceed  10 
pA. 

A.4.1.2  In-Source,  In-Situ,  and  Post-Irradiation  Measurements.  The  “post” 
irradiation  characterization  may  be  performed  in-source,  in-situ,  or  at  a  test  setup 
away  from  the  radiation  source.  The  details  of  performing  the  irradiation  are  pro¬ 
vided  in  MIL-STD-883  test  method  1019.3  and  MIL-HDBK-279,  “Piecepart  Total 
Dose  Hardness  Assurance  Guideline.”  The  test  procedure  for  measuring  the  Ip(s) 


76 


versus  Vas  characteristic  is  the  same  for  post-irradiation  as  for  pre-irradiation  with 
the  following  cautions.  First,  the  range  of  the  Vqs  sweep  to  go  from  leakage  through 
inversion  will  change  for  two  reasons:  a)  hole  trapping  in  the  oxide  will  cause  the  entire 
curve  to  shift  toward  negative  values,  and  b)  interface  states  will  cause  a  stretchout  of 
the  subthreshold  region.  Therefore,  the  range  of  Vqs  values  will  have  to  be  adjusted  to 
assure  that  the  entire  Id(S)  versus  V^s  characteristic  is  measured.  The  second  caution 
concerns  the  measurement  duty  cycle  for  in-source  measurements.  During  the  Id(s) 
versus  Vqs  measurement  cycle,  the  gate  voltage  is  usually  being  swept  from  negative 
to  positive  values  (for  n-channel  device  or  positive  to  negative  for  a  p-channel).  Since 
the  ionization  response  is  dependent  on  the  oxide  electric  field  magnitude  and  direc¬ 
tion  during  irradiation,  the  time  period  of  the  measurement  must  be  minimized.  The 
restrictions  on  measurement  time  will  depend  on  dose  rate,  the  dose  values  selected 
for  measurement,  and  the  total  number  of  test  samples  being  measured.  The  duty 
cycle  should  be  maintained  at  <1%. 

The  third  caution  concerns  temperature.  If  the  test  sample  temperature 
is  not  controlled  by  a  heating  or  cooling  stage,  the  temperature  will  increase  as  a 
result  of  the  irradiation.  If  the  measurements  are  in-situ,  either  the  dose  rate  must 
be  kept  low  enough  to  cause  an  insignificant  temperature  increase  in  the  MUT,  or  the 
radiation  source  and  MUT  must  be  cooled  with  forced  air  or  nitrogen.  For  in-situ 
measurements,  either  the  radiation  source  and  MUT  should  be  cooled,  or  the  MUT 
should  be  allowed  to  cool  prior  to  characterization  when  the  source  is  removed. 

With  each  set  of  measurements,  the  individual  /d(s)>Vgs  data  points  shall 
be  recorded  along  with  the  part  identification  number,  temperature  dose  rate,  dose 
and  time  between  °nd  of  irradiation  and  initiation  of  test. 


A.4.2  Data  Analysis. 

A.4.2.1  Determination  of  VTIj  and  A Vth>  V Tn  shall  be  determined  using  an 
extrapolation  of  the  drain  (or  source)  current  versus  the  gate-source  voltage.  If  mea¬ 
surements  were  made  using  Vps  =  100  —  150mV  (linear  operation),  a  plot  of  Id(s) 
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versus  Vqs  shall  be  made.  The  straight  line  portion  of  Id(s)  shall  be  extrapolated  to 
Id{s)  =  0.  Vth  is  the  Vqs  intercept  of  this  extrapolation.  If  the  value  of  Vds  was 
chosen  for  saturated  operation  then  a  plot  of  ^JId(s)  versus  F35  shall  be  made.  Again, 
the  straight  line  portion  of  \JId(s)  shall  be  extrapolated  Id(s)  =  0.  VTjj  is  the  Vgs 
intercept  with  this  extrapolation. 

The  value  of  AVjjy  is  found  from  the  following  relation, 


A  Vth  =  Vthi  -  Vthi 

where  Vth*i  is  found  from  the  post-irradiation  Id  versus  Vqs  characteristic,  and  VTHI 
is  found  from  the  pre-irradiation  characteristic. 


A.4.2.2  Calculation  of  I  mg  and  Iiw  The  values  of  mid^ap  current  I  mg  and 
inversion  current  Ijnv  are  found  from  the  following  equations. 
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N  =  silicon  doping  density  at  Si02  -  Si  interface 
n,-  =  intrinsic  carrier  concentration 
e t  =  silicon  dielectric  constant 
Cox  =  oxide  capacitance  per  unit  area 
C  =  ia*. 

•ox 

eox  =  oxide  (insulator)  dielectric  constant 

tox  =  oxide  thickness 

Z  =  channel  width 

Leff  =  effective  channel  length 

H  =  channel  mobility 

In  order  to  calculate  I  mg  and  IinVi  it  is  necessary  to  determine  the  quanti¬ 
ties  N,  tox  and  GM/VDS.  Table  1  is  a  list  of  methods  for  determining  these  parameters. 

Tablo  1.  Parameter  required  for  calculation. 

Parameter  Method  of  Determination 

a)  Vth  versus  substrate  bias. 

b)  Process  code  calculation  (e.g.,  SUPREM). 

c)  Spreading  resistance  profiles. 

a)  Ellipsometry. 

b)  C-V  Measurements. 

c)  Interferometry. 

a)  Plot  of  A/xj/AVg-s* 


N 


2m. 

VDS 


The  quantities  N  and  tox  are  normally  determined  by  the  manufacturer  and 


should  be  requested  for  each  sample  type  and  lot.  The  value  of  t,M ; 


(Vos  obtamc . 


from  the  inversion  data  of  Id{s)  versus  Vqs •  A  plot  of  AId(s)/ ^Vgs  versus  Vos  is 
made  such  as  the  one  shown  in  Figure  25.  The  value  of  Gm  used  in  the  calculation 
should  be  the  maximum  value  of  Gm  which  occurs  near  the  threshold  voltage. 
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GRAPHICS  PLOT  ****** 


Figure  25.  Plot  of  G\i  versus  gate  voltage  in  the  linear  region  for  drain 
voltage  of  0.1. 


Once  the  values  of  N,  tox,  and  Gm/Vds  are  obtained  for  each  MUT,  the 
values  of  Img  and  Iinv  can  be  calculated.  The  pre-irradiation  value  of  I  mg  and  Iinv 
should  be  used  for  all  calculations  involving  these  terms.  The  rationale  for  this  is 
discussed  in  the  guidelines. 


A.4.3  Determination  of  Vinv  and  AV/jvv* 

The  inversion  voltage,  Vinv,  is  defined  as  the  value  of  gate  to  source  voltage 
corresponding  to  the  calculated  value  of  Iinv-  The  value  of  Vjnv  is  found  from  a 
plot  of  log  Id(s)  versus  Vqs  by  interpolating  between  the  two  closest  measurements 
on  either  side  of  Iinv-  A  logarithmic  interpolation  is  preferred.  The  value  of  A Vinv 
is  given  by  the  relation 

A  Vinv  =  Vinv^  ~  Vinvi 

where  VINVl  is  the  post-irradiation  value  of  Vinv  and  Vinvi  is  the  pre-irradiation 
value. 


A. 4. 4  Determination  of  A Vmg  and  AS. 

In  order  to  apply  the  subthreshold  charge  separation  technique,  it  is  nec¬ 
essary  to  determine  AVWG  and/or  AS  from  the  pre-  and  post-irradiation  subthresh¬ 
old  characteristics.  A  Vmg  is  the  difference  between  the  post-irradiation  value  of  the 
midgap  voltage,  Vmg'i,  and  the  pre-irradiation  value,  Vmgi-  The  values  of  Vmgi  and 
Vmgi  are  found  from  a  plot  of  log  Ip(s)  versus  Vqs  for  the  pre-  and  post-irradiation 
data  set  of  Id(s),Vgs  points.  A  typical  set  of  pre-  and  post-irradiation  curves  are 
shown  in  Figure  26.  The  straight  line  portion  of  the  subthreshold  log  Id(S)  versus  Vqs 
curve  is  extrapolated  to  the  calculated  value  of  Img-  The  value  of  VG5  at  Img  is  the 
value  of  Vmg-  For  each  post-irradiation  data  set 


A  Vmg  =  Vmg-i  -  Vmgi- 
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The  subthreshold  swing,  A  S',  is  equal  to  the  post-irradiation  value  of  the 
subthreshold  swing,  S7,  minus  the  pre-irradiation  value  Si.  The  values  of  Si  and 
are  found  from  the  same  straight  lines  used  to  extrapolated  to  the  values  of  Vmgi  and 
VxfG',  as  shown  in  Figure  26.  The  values  S'/  and  S7  are  measured  in  volts  per  decade 
of  current. 


A.4.5  Separation  of  A Vth  or  AVinv  into  AV0T  and  A V/r. 

The  charge  separation  technique  can  be  applied  to  either  the  threshold  volt¬ 
age  shift,  A  Vth,  or  the  inversion  voltage  shift,  AV1NV.  A  Vth  >s  often  preferred  since 
it  is  the  parameter  typically  specified  and  measured  as  a  performance  parameter. 
However,  A ViNV  has  greater  physical  significance  since  it  is  well  defined  in  terms  of 
the  surface  potential.  Although  either  parameter  may  be  used  for  charge  separation, 
A VTh  is  usually  preferred  for  engineering  applications. 


A. 4. 5.1  Method  A:  Midgap  Voltage  Techniques.  In  the  midgap  voltage  tech¬ 
nique  the  value  of  AVot>  the  voltage  component  due  to  oxide  trapped  charge,  is  found 
from  the  value  of  midgap  voltage  shift, 


AVqt  =  AVmg  =  VgM'i  -  Vmgi 


The  value  of  A  V/y,  the  voltage  component  due  to  interface  states,  is  then  equal  to  the 
difference  between  AVth{AVihv)  and  AVor>  be., 


A  Vit  =  A  Vth  (AVihv)  -  AVqt 


A.4.5. 2  Method  B:  Subthreshold  Swing  Technique.  In  the  subthreshold  swing 
technique,  the  value  of  AVn  is  found  from  the  value  of  AS  according  to  the  following 
relation 


AV it  = 


(  q<^B  \ 

\kTlnlo) 


AS 
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where  AS  —  S^-  Si  and  <f>B  =  kT/qln(N/rii).  The  value  of  A Vot  is  then  found  from 
AVth(AVinv)  according  to  the  relation 


AVqt  =  A Vth  (AV/^rv)  _  AV/j* 


A. 5  REPORTING. 

Charge  separation  reporting  shall,  at  a  minimum,  include  the  following  in¬ 
formation:  a)  part  identification;  b)  irradiation  data  including  dose  rate,  dose  levels, 
irradiation  bias,  anneal  bias,  time  between  irradiation  and  test,  test  type  (in-source, 
in-situ  or  pre/post)  and  test  temperature;  c)  values  of  N  and  tox  used  in  the  calcula¬ 
tions;  d)  whether  Id  or  Is  was  used;  e)  VDs\  f)  whether  AVTH  or  AVinv  was  used; 
g)  data  separation  method  (A  or  B);  and  h)  calculated  values  of  AVOT  and  AV;T  at 
each  dose. 


A. 6  SUMMARY. 

The  following  details  shall  be  specified  in  the  applicable  procurement  docu¬ 
ment. 

1.  Part  types  and  quantities  to  be  tested. 

2.  Test  temperature. 

3.  Test  instruments  if  other  than  specified  in  A. 2.1. 

4.  Irradiation  test  methods,  dose  rate,  dose  levels,  irradiation  and  anneal  biases, 
and  measurement  times. 

5.  VDs  for  linear  or  saturated  measurement. 

6.  Charge  separation  method  A  or  B. 

7.  Whether  A  Vth  or  AVinv  is  used  for  charge  separation. 
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